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EXECUTIVE  SUMMARY 


This  document  presents  an  annual  report  to  the  Office  of  Naval  Research  for  a  research  program 
witirterf  "A  Theoretical  Search  For  Supervelodty  Semiconductors".  This  program  has  been  funded  by 
ONR  since  1974  in  the  Department  of  Electrical  and  Computer  Engineering  at  N.C.  State  University. 
The  research  has  resulted  in  more  than  80  refereed  publications  and  numerous  conference  presentations 
from  its  inception.  Major  contributions  to  die  field  of  hot  electron  transport  and  semiconductor  device 
modeling  have  been  achieved,  new  computational  methods  have  been  developed  (e.g.  path  integral 
Monte  Carlo  techniques),  and  the  work  has  helped  atimulatc  commercial  ventures  in  die  applications  of 
quaternary  semiconductor  materials  to  electronic  and  optical  devices.  In  addition,  there  have  been 

twenty  eight  Ph.D.  and  M.S.  students  who  have  received  degrees  at  N.C.  State  University  with 

» 

research  support  from  this  contract  Three  visiting  faculty  members  horn  Japan  came  to  the  University*,** 
to  work  with  die  faculty  investigators  supported  under  this  ONR  contract  during  die  1979-1983  time 
period.  A  visiting  professor  from  die  French  CNRS  Microstructures  and  Microelectronics  Laboratory 
in  Bagneux  (near  Paris)  spent  a  sabbatical  year  at  N.C.  State  during  1988-89,  and  he  devoted  full-time 
working  on  this  program  at  no  cost  to  ONR.  During  die  current  funding  period,  a  visiting  scholar 
from  China  is  a  member  of  our  research  group  working  on  projects  which  directly  impact  this  ONR 
program.  This  researcher  will  be  appointed  as  a  Visiting  Assistant  Professor  in  August,  1992,  and  will 
spend  approximately  30%  of  his  time  working  on  this  project  His  work  will  be  devoted  to  advanced 
device  simulation  and  physics.  In  particular,  he  will  develop  new  concepts  for  die  real  space  transfer 
logic  transistor  and  other  charge  injection  transistor  structures.  Finally,  we  have  hired  a  new  post¬ 
doctoral  research  associate  from  the  University  of  Michigan  (Prof.  J.  Singh’s  group)  and  he  will  join 
us  to  work  on  transport  theory  in  various  quantum-based  devices. 

This  initial  phases  of  this  work  centered  around  die  development  of  Monte  Carlo  simulation  tech¬ 
niques  which  allow  the  study  of  detailed  physics  of  hot  electron  transport  in  a  variety  of  compound 


semiconductor  materials.  The  original  emphases  were  concerned  with  electronic  materials  phenomena. 
Later  work  considered  the  utilization  of  these  materials  in  realistic  device  structures  where  physical 
boundary  conditions  must  be  imposed  on  the  carrier  transport  More  recently,  the  work  has  focused  on 
the  domain  of  ultra-small  materials  and  device  phenomena  where  microscopic  non-local  transient 
effects  such  as  velocity  overshoot  ballistic  and  neariy-ballistic  transport  and  quantum  transport 
become  important  or  dominant  During  the  past  five  years  we  have  researched  die  applications  of  the 
Feynman  "integral  over  paths”  approach  to  quantum  transport  and  identified  numerical  limitations  to  its 
practical  application.  Also  we  have  emphasized  the  study  of  hot  electron  effects  in  new  device  struc¬ 
tures,  such  as  the  hot  electron  spectrometer,  heterojunction  bipolar  transistor,  small  dimension  metal- 
semiconductor-metal  photodetectors,  delta-doped  high  electron  mobility  transistors,  and  real  space 
transfer  logic  transistors.  During  this  time,  we  have  been  exploring  some  new  approaches  to  device 
modeling  which  combine  the  Monte  Carlo  method  with  the  method  of  moments  of  die  Boltzmann  tran¬ 
sport  equation  (hydrodynamic  transport  model)  for  studying  specific  device  structures,  such  as  small- 
dimension  n+-n-n+  majority  carrier  devices  and  the  high  electron  mobility  transistor.  In  addition,  we 
have  incorporated  quantum  correction  terms  into  the  hydrodynamic  model  and  applied  this  model  to 
resonant  tunneling  structures.  Recently,  we  have  applied  a  new  ansaiz  distribution  function  as  a  consti¬ 
tutive  relation  to  dose  the  moment  equations  in  die  hydrodynamic  transport  model.  Initial  results  of 
this  approach  have  been  physically  satisfying  and  computationally  promising.  A  comprehensive  formu¬ 
lation  of  this  new  modd  has  been  published  and  has  attracted  a  lot  of  attention. 

During  the  next  year,  new  research  in  four  general  areas  will  advance  our  basic  research  direc¬ 
tions  encompassing  the  study  of  hot  electron  transport  in  materials  and  devices.  These  four  areas 
included  1)  the  study  of  quantum  transport  in  mesoscopic  structures,  nanostructures,  and  related  ultra¬ 
small  heterobarrier  device  structures  using  several  physical  models  and  computational  approaches,  2) 
new  approaches  to  the  merging  of  Monte  Carlo  methods  with  moment  equation  methods  with  resulting 
improvements  to  the  hydrodynamic  transport  model,  3)  Monte  Carlo  simulations  in  order  to  study 


pseudomorphic  devices,  real  space  transfer  stmctures,  and  submicron  MOSFETs,  and  4)  search  for  new 
quantum-based  device  concepts  for  ulna-fast,  ultra-dense  applications.  The  possibility  of  expanding 
our  experience  to  optical  and  optoelectronic  device  modeling  will  be  seriously  explored  as  well,  as  we 
define  the  role  of  our  new  personnel  which  will  be  added  to  our  group  in  Fall  1992. 
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1.0  INTRODUCTION 

In  October,  1974,  the  Office  of  Naval  Research  initiated  sponsorship  of  a  basic  research  program 
in  die  Department  of  Electrical  and  Computer  Engineering  at  North  Carolina  State  University.  The 
general  goal  of  this  research  program  has  been  and  continues  to  be  die  investigation  of  high  speed  car¬ 
rier  transport  in  DI-V  compound  semiconductors,  m-V  alloy  materials,  other  advanced  electronic  and 
optical  materials,  and  novel  device  structures  which  utilize  these  materials.  Four  faculty  members  at 
N.  C  State  have  been  primarily  involved  and  supported  by  this  project  Currently,  Profs.  M.  A. 
Littlejohn  and  K.  W.  Kim  serve  as  the  co-principal  investigators  on  tins  project  This  research  pro¬ 
gram  has  made  significant  contributions  to  the  understanding  and  knowledge  base  of  hot  electron  tran¬ 
sport  in  materials  and  devices.  It  has  provided  scientific  guidance  to  the  U.  S.  Navy  in  the  formulation 
of  a  part  of  its  basic  and  applied  research  program.  Numerical  concepts  developed  under  this  project 
have  been  transferred  to  other  Universities,  jnrfnHtwg  the  University  of  Illinois,  and  we  continue  these 
inter-institutional  collaborations.  During  the  past  year,  we  have  maintained  dose  contact  with  die 
research  program  of  Prof.  Hess  at  die  University  of  Illinois.  Also,  collaborations  with  Dr.  H.  L.  Gni- 
bin  of  Scientific  Research  Associates,  Inc.  and  Dr.  G.  J.  Iafrate  and  Dr.  M.  A.  Strosdo  of  the  U.  S. 
Army  Research  Office  in  the  area  of  quantum  transport  in  semiconductor  devices  have  been 
strengthened.  Dr.  Iafrate,  Dr.  Strosdo,  and  Dr.  Grubin  are  Adjunct  Professors  in  the  ECE  Department 
at  North  Carolina  State  University  and  we  share  graduate  students  and  post-doctoral  research  associates 
on  joint  projects.  In  addition  to  those  directions,  our  research  results  have  helped  stimulate  commercial 
ventures,  particularly  in  die  development  of  GalnAsP-based  materials  and  devices.  To-date,  this  pro¬ 
gram  has  resulted  in  79  refereed  publications  in  the  literature,  1  additional  manuscript  is  currently  in 
press,  and  6  manuscripts  have  been  submitted  or  are  in  preparation  for  submission  to  die  technical 
literature.  A  listing  of  these  publications  is  given  in  Appendix  A.  In  addition,  numerous  invited  talks 
and  presentations  have  been  given  at  conferences  and  workshops  throughout  the  United  States  and  in 
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other  countries.  The  program  has  contributed  significantly  to  the  educational  program  at  N.  C.  State 
University  with  more  than  thirty  Ph.D.,  MS.  and  undergraduate  students  having  received  support  under 
this  ONR  contract  Currently,  three  PhJX  students  are  working  toward  their  degrees  on  this  project 
Two  of  these  Ph.D.  students  are  U.S.  Citizens.  A  visiting  scholar  from  China,  who  will  be  appointed 
as  a  Visiting  Assistant  Professor  in  August  1992,  is  also  working  on  research  related  to  this  program. 
We  have  hired  another  post-doctoral  research  associate  who  will  extend  our  capabilities  into  the  area  of 
band  structure  effects,  quantum  transport,  and  many-body  physics. 

This  ONR  program  has  been  efficient  and  productive.  The  high  quality  of  this  research  will  con¬ 
tinue  to  be  maintained  in  die  future.  This  report  summarizes  the  progress  and  accomplishments  made 
during  the  past  contract  period. 
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2J>  RESEARCH  RESULTS 


2.1  Background 

Since  its  inception,  semiconductor  technology  has  been  stimulated  by  requirements  for  electronic 
systems  with  ever-increasing  capabilities  to  process  information  faster,  more  functionally  and  more 
efficiently.  These  requirements  have  motivated  foe  scaling  down  of  integrated  circuit  (IQ  device 
dimensions  into  foe  submicron  (less  than  ten  thousand  angstroms)  and  ultrasubmicron  (less  than  one 
thousand  angstroms)  regions.  Today,  we  have  entered  an  era  where  nanostructure  physics  and  fabrica¬ 
tion  motivate  our  research  efforts  in  semiconductor  electronics  [1].  As  fabrication  technology  has 
allowed  such  devices  to  be  realized,  many  new  and  fundamental  questions  have  emerged  concerning 
the  underlying  physics  of  small  (atomic  level)  dtm«Mrinm  in  semiconductor  devices.  Important  issues 
now  under  consideration  for  ultrasubmicron  devices  include  nonequilibrium  transport  dealing  with  such 
topics  as  quasi-ballistic  transport,  overshoot  phenomena  and  quantum  transport.  A  great  deal  of  pro¬ 
gress  has  been  achieved  in  our  understanding  of  these  important  device  effects,  although  major  work 
remains  to  be  done  as  our  ability  to  fabricate  very  small  electronic  device  structures  continues  id 
expand  and  mature  [2,3]. 

The  ability  to  fabricate  small  devices  has  been  continually  refined  over  foe  last  ten  years  through 
impressive  improvements  in  materials  growth  technologies.  Molecular  beam  epitaxy  (MBE),  metalor- 
ganic  chemical  vapor  deposition  (OMCVD),  and  atomic  layer  epitaxy  (ALE)  have  provided  foe  ability 
to  fabricate  a  wide  variety  of  materials  and  heterostructure  combinations  with  near  perfect  interfaces, 
doping  control  and  compositional  uniformity  with  atomic  level  dimensions.  The  development  of  ALE 
may  very  well  prove  to  be  the  ultimate  growth  technology  since  it  allows  foe  deposition  of  one  mono¬ 
layer  of  device  quality  materials  through  a  controllable,  self-limiting  mechanism,  and  is  especially  use¬ 
ful  for  foe  deposition  of  heterojunctions  [4].  The  ability  to  grow  layers  with  dimensions  of  a  few 
angstroms  opens  the  domain  of  quantum  transport  to  experimental  study  and  verification.  Thus,  topics 
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resulting  from  size  quantization  in  condensed  matter  must  be  investigated  from  theoretical  viewpoints 
with  tools  which  are  either  partially  developed  or  through  the  development  of  new  tools  which  are  not 
now  available.  Quantization  effects  arising  from  geometrical  size  constraints,  proximity  effects  result¬ 
ing  from  closely  packed  arrays  of  devices,  and  general  solid-state  considerations  not  heretofore  con¬ 
sidered  questionable  (effective  mass  approximation,  the  role  of  contacts  and  die  like)  must  be 
addressed  from  a  fundamental  point  of  view.  Moreover,  from  a  device  physics  point-of-view,  it  is 
desirable  to  have  a  microscopic  description  of  the  physics  of  smalt  dimensions  which  is  amwiahi*  to 
phenomenological  treatment,  so  that  its  properties  can  be  meaningfully  incorporated  into  futuristic  dev¬ 
ice  concepts  and  simulations. 

Theoretical  methods  to  address  carrier  transport  have  also  progressed  rapidly  over  the  last  ten- 
fifteen  years,  in  a  manner  similar  to  research  in  semiconductor  thin  film  epitaxial  growth  technology. 
In  fact,  this  is  a  natural  progression  in  many  ways  and  is  to  be  expected.  The  progress  achieved  in 
materials  growth  of  structures  with  quantum  dimensions  dictates  that  new  approaches  be  developed  and 
refined  to  study  quantum  transport  phenomena  and  the  physics  of  small  dimensions.  However,  a 
significant  change  in  direction  is  now  warranted.  Past  transport  theory  and  device  modeling 
approaches  have  relied  on  particle  or  quasi-particle  approaches  where  the  electrons  are  treated  as  rigid 
mobile  entities  which  undergo  interactions  with  the  transport  medium.  The  treatment  of  the  interaction 
often  involves  wave  concepts.  However,  the  model  is  basically  a  particle  model.  In  the  current  regime 
of  quantum  transport,  we  may  no  longer  be  able  to  consider  the  carriers  as  particles.  It  is  quite  likely 
that  their  physical  behavior  will  be  governed  either  partially  or  completely  by  wave  phenomena. 

Quasi-particle  methods  are  an  attempt  to  model  structures  with  quantum  dimensions  which  retain 
as  much  of  the  classical  formalism  as  possible  in  order  to  be  able  to  express  results  in  terms  of  param¬ 
eters  which  are  of  tire  greatest  experimental  interest,  such  as  carrier  velocity  and  diffusion  constant. 
This  is  a  flexible  approach.  However,  much  care  is  required  to  ensure  that  all  important  effects  are 
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property  included  because  of  die  approximations  involved  in  die  formulations.  On  the  other  hand, 
more  fundamental  quantum  approaches,  such  as  operator-eigenfunction  methods,  adhere  closely  to  the 
actual  quantum  states  present  in  the  device  structure  when  scattering  is  not  included.  Scattering 
processes  (dissipation)  can  be  added  by  using  perturbation  theory  from  quantum  mechanics.  These 
techniques  can  obtain  the  greatest  sensitivity  to  the  resulting  carrier  confinement  and  the  lattice  poten¬ 
tials.  However,  they  are  relatively  inflexible  in  studying  non-linear  dynamical  properties  in  the  pres¬ 
ence  of  strong  dissipation,  such  as  is  present  in  the  electron-phonon  interaction  at  high  electric  fields. 

Another  approach  to  quantum  transport  relies  cm  the  "integral  over  paths"  method,  originally  pro¬ 
posed  by  Dirac  and  formulated  by  Feynman  [5].  Practical  path  integral  methods  for  die  study  of  small 
devices  rely  on  an  influence  functional  technique  in  which  the  source  of  the  dissipation  has  been 
integrated  over  all  phonon  modes.  This  results  in  a  model  influence  functional  where  the  phonon- 
scattering  dissipation  can  be  represented  as  an  interaction  with  a  collection  of  harmonic  oscillator 
modes  in  which  die  translational  invariance  of  die  carriers  is  preserved.  The  resulting  model  includes 
constant  or  oscillatory  electric  and  magnetic  fields,  carrier  screening,  scattering  and  dissipation,  carrier 
confinement,  background  temperature  and  initial  conditions  can  be  dealt  with  as  readily  as  for  a  free 
particle  [6].  The  path  integral  method  can  be  compared  and  contrasted  to  other  methods  now  being 
studied  and  supported  by  ONR,  such  as  the  application  of  the  density  matrix  formalism  and  Wigner 
distribution  function  approaches  including  the  use  of  moment  equations  [7].  Based  on  our  progress  in 
developing  the  path  integral  approach  during  the  past  five  years,  we  believe  that  path  integral  methods 
will  play  an  important  role  in  developing  an  understanding  of  quantum  transport  in  reduced  geometries 
(such  as  mesoscopic  devices  {8]).  However,  because  of  its  intensive  computation  requirements,  the 
path  integral  method  is  very  limited  in  its  application  to  practical  devices.  Currently,  we  are  aware  of 
only  one  other  U.S.  university  program  directed  toward  applying  the  path  integral  method  to  semicon¬ 
ductor  devices.  This  is  the  program  at  the  University  of  Illinois  under  the  direction  of  Prof.  Kart  Hess. 
Our  efforts  in  die  development  of  this  method  are  complimentary  to  the  me  at  Illinois,  and  we  have 
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collaborated  with  die  Illinois  personnel  through  the  use  of  the  NSF  Supercomputing  Center  and  the 
National  Center  for  Computational  Electronics  located  there.  One  of  our  former  students,  whose  Ph.D. 
dissertation  on  path  integral  methods  was  supported  by  this  ONR  contract,  joined  Prof  Hess'  group  in 
January.  1990  as  a  Visiting  Research  Assistant  Professor.  We  {dan  to  maintain  contact  with  Prof. 
Hess'  group  as  we  evaluate  our  path  integral  research  directions.  We  win  continue  this  work  if  we  can 
hire  an  appropriate  Ph.D.  candidate. 

The  work  at  N.  C.  State  University  supported  by  the  ONR  under  the  current  contract  has  pro¬ 
gressed  over  the  past  fifteen  years  from  the  realm  of  particle  models  to  quasi-particle  models  to  quan¬ 
tum  transport  models.  We  have  relied  primarily  an  the  Monte  Carlo  method  to  study  and  solve  tran¬ 
sport  problems  in  III-V  compound  semiconductors.  Our  techniques  for  modeling  materials  physics  and 
device  phenomena  extend  to  device  dimensions  around  one  thousand  angstroms  [10.11],  and  we  are 
confident  in  these  models  for  predicting  steady  state  and  transient  device  effects  down  to  these  dimen¬ 
sions.  However,  in  order  to  remain  in  the  forefront  of  transport  physics  and  device  research,  we  must 
continue  and  increase  our  progress  into  the  realm  of  dimensions  where  physical  effects  can  be  studied 
in  device  regions  with  dimensions  less  than  one  thousand  angstroms.  We  will  continue  to  explore  the 
Monte  Carlo  transport  method  for  study  of  novel  device  structures  where  theoretical  underpinning  is 
required.  One  novel  class  of  structures  which  will  be  studied  in  detail  by  the  Monte  Carlo  method  is 
based  on  the  rapidly-emerging  pseudomorphic  or  strained-layer  devices.  These  structures  permit 
extended  compositional  ranges  and,  thus,  have  a  number  of  potential  advantages  such  as  higher  tran¬ 
sconductance  and  channel  carrier  density.  Recent  studies  on  the  strain-induced  piezoelectric  fields  have 
opened  yet  mother  possibility  in  realizing  ultra-fast  switching  devices,  resulting  in  an  increased  impor¬ 
tance  on  the  study  of  carrier  transport  in  these  devices.  Other  novel  devices  which  will  be  analyzed 
include  delta-doped  field-effect  translators  and  real-space  transfer  devices,  based  on  studies  which  have 
begun  during  the  past  two  years.  In  addition,  we  will  incorporate  Monte  Carlo  methods  into  the 
quasi-particle  approach  based  on  moments  of  the  Boltzmann  transport  equation  with  quantum 
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mechanical  corrections.  Preliminary  results  show  the  utility  of  this  approach  for  modeling  microwave 
and  millimeter  wave  devices  which  are  important  to  the  DoD’s  MMIC  program. 

Finally,  we  want  to  re-iterate  the  important  interactions  which  have  developed  between  our 
research  group  and  other  researchers  during  the  last  year.  These  interactions  increase  the  impact  of 
this  pmginn  on  the  field  of  semiconductor  device  physics  and  allow  our  research  efforts  to  be  far  more 
productive  through  increased  intellectual  efforts  and  enhanced  facilities  and  resources.  During  the 
1991-1992  funding  period,  we  have  maintained  collaborations  with  Prof.  K.  Hess  at  the  University  of 
Illinois,  Dr.  H.  Grubin  of  Scientific  Research  Associates,  and  Drs.  M.  Stroscio  and  G.  Iafrate  of  the 
U.  S.  Army  Research  Office.  These  research  collaborations  have  resulted  in  several  joint  publications 
based  on  mutual  research  interests,  expertise  and  capabilities.  There  have  been  several  visits  between 
these  laboratories  and  die  logistics  for  increased  collaborations  are  excellent 

Thus,  we  believe  that  we  have  a  well-rounded  approach  for  this  ONR  program  which  will  lead  to 
new  innovations,  and  practical  developments  in  semiconductor  device  physics,  simulation  and  model¬ 
ing. 

12  Summary  of  Research  Results 

This  section  provides  die  current  status  of  our  research  efforts  in  specific  topical  areas  we  have 
worked  during  the  past  year. 

2.2, 1  Hydrodynamic  hot-electron  transport  model 

The  primary  goal  of  this  research  is  to  develop  and  implement  an  advanced  hot-electron  transport 
model  suitable  for  accurate  and  ef  lciem  study  of  electron  dynamics  in  both  submicron  and  ultrasmall 
semiconductor  device  structures.  The  model  is  required  to  accurately  reflect  hot-electron  effects  such 
as  velocity  overshoot  and  nonstationary  effects  due  to  spatial  nonuniformities  in  average  election  con¬ 
centration,  avenge  electron  velocity  aid  average  election  energy.  Also,  die  model  is  required  to  accu- 
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ratcly  describe  die  streaming  motion  and  dissipation  of  electrons  in  multiple  ncmparabolic  conduction 
bands  (as  are  present  in  especially  all  m-V  compound  semiconductors).  In  addition,  an  election  tran¬ 
sport  model  is  desired  which  can  be  extended  to  the  study  of  quantum  systems  (i.e.,  feature  sizes  less 
than  1000  A)  where  wave-lQce  characteristics  (e.g.,  tunneling  and  interference  effects)  are  important. 

Over  the  past  decade,  much  effort  has  been  focused  on  the  development  of  semi -classical  electron 
transport  models  which  accurately  reflect  hot-electron  effects  in  submicron  semiconducting  device 
structures  [11].  The  incentive  is  that  nonequilibrium  and/or  nonstationary  effects  may  be  employed  to 
achieve  superior  electron  device  performance.  One  important  example  of  a  potentially  useful  hot- 
electron  effect  is  velocity  overshoot  [12].  In  velocity  overshoot,  the  difference  in  momentum  and 
energy  relaxation  times  can  result  in  an  initial  temporal  drift  velocity,  under  constant  electric  field  con¬ 
ditions,  which  overshoots  the  corresponding  time  steady-state  (stationary)  value.  This  dynamical 
response  in  electron  motion  is  due  to  the  difference  between  the  displacement  of  the  electron  distribu¬ 
tion  in  momentum  space  and  the  displacement  in  energy  space  (i.e.,  the  form  of  the  nonequilibrium 
distribution  function).  For  this  type  of  electron  motion,  (me  cannot  use  simple  models  based  upon  sta¬ 
tionary  relationships  (which  ignore  important  features  of  the  distribution)  to  describe  the  transport  phy¬ 
sics.  Furthermore,  deficiencies  in  these  near-equilibrium  models  may  also  fail  to  accurately  describe 
electron  transport  in  submicron  device  structures  even  in  time  steady-state.  This  can  occur  when  die 
electron  gas  is  far  from  equilibrium  and  large  spatial  gradients  in  the  electron  gas  parameters  (electron 
density,  average  electron  velocity,  average  electron  energy,  etc.)  exist 

Since  classical  (drift-diffusion)  models  assume  local  stationary  relationships  (average  velocity 
depends  on  the  local  electric  field)  and  are  based  upon  near  equilibrium  derivations  [12],  physically 
superior  models  must  be  used  to  fully  and  effectively  study  important  hot-electron  effects.  At  present 
the  problem  of  treating  transport  in  submicron  structures  is  typically  attacked  using  one  of  two  distinct 
methods.  The  first  and  most  accurate  approach  has  been  to  apply  Monte  Carlo  simulation  methods 
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[13].  The  Monte  Carlo  method  is  a  powerful  statistics-based  numerical  technique  for  investigating 
electron  transport  physics  in  semiconductors  and  semiconductor  devices  [13].  While  the  Monte  Carlo 
method  provides  a  relatively  simple  and  accurate  indirect  approach  to  determine  panicle  distribution 
functions  and  transport  parameters  resulting  from  the  Boltzmann  transport  equation  subject  to  complex 
«ft»w.ring  mechanisms,  thr  approach  does  have  important  disadvantages.  Specifically,  the  technique 
often  requires  extensive  computation  time  to  arrive  at  adequate  statistics  and  does  not  naturally  possess 
macroscopic  terms  for  physical  interpretation. 

The  alternative  (second)  approach  to  study  electron  device  physics  more  accurately  and  efficiently 
is  to  develop  macroscopic  models  which  predict  average  transport  quantities  as  opposed  to  the  full 
electron  distribution  function.  In  this  approach,  foe  description  of  electron  dynamics  is  reduced  to 
hydrodynamic-like  conservation  equations  for  electron  density,  average  electron  momentum,  average 
electron  energy,  etc  [14].  The  hydrodynamic  models  are  derived  horn  moments  of  foe  Boltzmann  tran¬ 
sport  equation  which  usually  indude  many  simplifying  assumptions  (i-e..  usually  neglect  conduction 
band  nonparabolicity,  treat  multiple  conduction  band  valley  effects  indirectly  and  directly  or  indirectly 
make  incorrect  assumptions  about  foe  form  of  foe  electron  distribution  function)  about  foe  electron 
transport  physics.  However,  if  advanced  hydrodynamic  models,  which  indude  hot  electron  effects 
such  as  multiple  nonparabolic  conduction  bands,  are  developed  and  combined  with  Monte  Carlo  tech¬ 
niques  this  is  an  excellent  overall  approach.  A  limited  number  of  Monte  Carlo  simulations  can  then  be 
utilized  to  verify  foe  accuracy  of  foe  simplified  models  and  to  supply  approximate  parametric  dissipa¬ 
tion  data  for  foe  particular  materials  and/or  devices  under  consideration.  These  types  of  advanced 
hydrodynamic  models,  together  with  the  Monte  Carlo  method,  represent  a  simplified  mathematical 
approach  for  accurate  and  efficient  study  of  realistic  transport  physics. 

Finally,  very  few  hydrodynamic  (macroscopic  or  average  ensemble-electron)  models  have  been 
implemented  to  accurately  represent  quantum  effects  (i.e.,  tunneling  and  interference)  present  in 
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ultrasmall  electron  structures.  One  method  to  accurately  model  structures  with  spatial  dimensions  on 
the  scale  of  the  de  Broglie  wavelength  (Le.,  in  the  quantum  regime)  is  to  utilize  quantum  hydrodynam¬ 
ics  [15,16].  In  this  method,  a  set  of  quantum  moment  equations  are  derived  by  taking  moments  [17] 
of  die  Wigner-Bohzmann  equation  and  closure  is  obtained  using  an  approximate  nonequflibrium  form 
of  die  zero-current  (equilibrium)  Wigner  distribution  function.  This  procedure  results  in  a  set  of 
quantum-corrected  hydrodynamic  transport  equations  which  have  been  demonstrated  to  predict  quantum 
phenomena  [15,16]. 

During  the  past  year,  advanced  forms  of  the  hydrodynamic  transport  model  (HTM)  were 
developed  and  utilized  to  study  both  submicron  and  ultrasmall  semiconducting  device  structures.  First 
an  advanced  semi-classical  hydrodynamic  transport  model  was  developed  and  implemented  to  accu¬ 
rately  predict  election  dynamics  in  multiple  nonparabolic  energy  bands.  In  addition,  a  new  quantum 
mechanical  formulation  of  the  first  three  moment  equations  (Le.,  quantum  hydrodynamic  equations  for 
electron  density,  average  velocity  and  average  energy)  was  investigated  and  implemented  to  predict 
non-classical  tunneling  in  a  double-barrier  heterostructure.  The  development  of  the  advanced  model(s) 
was  achieved  through  several  separate  theoretical  investigations. 

First,  to  determine  a  method  to  include  the  effects  of  conduction  band  nonparaboiicity  into  the 
macroscopic  streaming  (collision  independent)  of  the  semi -classical  HTM.  the  fundamental  physics  of 
electron  transport  were  investigated.  An  alternate  formalism  was  developed  for  deriving  a  new  nonpar¬ 
abolic  hydrodynamic  transport  model  adequate  for  modeling  hot  electron  transport  in  submicron  sem¬ 
iconductor  devices.  Specifically,  the  model  equations  were  developed  by  applying  a  unique  set  of 
moment  operators  ($o(k)=l;  a  constant,  0|(k)su(k);  the  nonparabolic  velocity, 
(m(k)/2]u(k)  u(k>»E(k);  an  approximation  to  the  particle  energy)  to  the  collisionless  Boltzmann  tran¬ 
sport  equation.  The  operators,  used  in  this  analysis,  represent  a  significant  improvement  in  die 
development  of  the  hydrodynamic  transport  model  since  they  lead  to  a  more  compact  mathematical 
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form.  This  resulting  form  can  be  manipulated  more  easily  and  reveals  dearly  the  allowable  simplify¬ 
ing  approximations.  An  intuition-based  distribution  function  was  proposed  as  a  constitutive  (essential) 
relation  to  dose  the  moment  equations.  This  ansatz  nonparabolic  distribution  function,  which  is 
derived  in  vdority  space  for  mathematical  convenience  and  approximates  some  of  the  characteristics  of 
die  true  nonparabolic  distribution  function,  was  then  analyzed  and  used  to  devdop  a  set  of  transport 
parameters.  After  development,  the  results  from  two  independent  Monte  Carlo  transport  models  were 
used  to  evaluate  the  nonparabolic  transport  parameters.  These  models  were  used  to  investigate  electron 
transport  in  both  stationary  (uniform  dectric  field)  and  nonstatkmary  (ballistic  diode)  environments. 
The  results  of  these  studies  show  that  there  is  excellent  agreement  between  tins  new  nonparabolic  tran¬ 
sport  model  and  the  Monte  Carlo  calculations.  Hence,  this  new  nonparabolic  transport  modd  presents 
a  viable  alternative  to  studying  realistic  electronic  device  structures  operating  under  high  bias  condi¬ 
tions. 

In  the  previous  analysis,  specific  macroscopic  collision  terms  are  omitted.  This  is  done  to  permit 
a  concise  focus  on  the  streaming  (collision  independent)  terms  that  evolve  from  the  application  of  non¬ 
parabolic  conduction  bands.  The  collision  terms,  which  are  almost  always  treated  phenomenologically, 
has  been  devdoped  separatdy  for  a  nonparabolic  multi-valley  system  (GaAs).  Specifically,  multi¬ 
valley  dissipation  is  included  in  the  modd  using  an  ensemble  relaxation  time  approximation.  The 
numerous  intravalley  and  intervalley  relaxation  times  were  calculated  using  tire  Monte  Carlo  Method. 

In  order  to  extend  the  modeling  approach  to  the  quantum  device  regime,  a  technique  was  investi¬ 
gated  to  devdop  quantum  corrections  for  the  hydrodynamic  equations.  The  set  of  quantum-corrected 
transport  equations  which  we  have  developed  snd  investigated  were  based  upon  quantum  moment 
equations  of  Grubin  and  Kreskovsky  [15].  Grubin  and  Kreskovsky  previously  constructed  a  set  of 
quantum  balance  (hydrodynamic)  equations  applicable  for  electron  transport  in  mesoscopic  structures. 
Their  model  was  derived  from  tire  general  moment  equations  of  Strosdo  [17]  using  the  displaced  none- 
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quflibrium  Wigner  distribution  function  of  Anoona  and  Iaftate  [18].  These  resulting  quantum  hydro¬ 
dynamic  equations  incorporate  quantum  effects  through  the  introduction  of  terms  (corrections)  with 
density  gradient  dependence. 

Both  the  advanced  semi-classical  and  the  quantum-corrected  hydrodynamic  transport  models 
(HTMs)  were  used  to  study  electron  transport  in  novel  semiconductor  device  structures.  A  general 
numerical  algorithm  was  implemented  to  solve  the  resulting  system  of  nonlinear  differential  equations 
[from  both  physical  models)]  for  arbitrary  one-space-dimensional  n-type  structures.  The  numerical 
algorithm  ntfflM*  a  local  simultaneous-coupled  (Newton)  iteration  approach  in  tandem  with  a  global 
physics-based  continuation.  Simulation  results  were  generated  for  a  single-valley  (T)  HTM,  with  quan¬ 
tum  corrections,  applied  to  a  double-barrier  heterostructure  with  A^Oa^As  barriers.  This  study  was 
die  first  to  use  the  three  quantum  hydrodynamic  (QHD)  equations  to  study  transport  in  heterostructure 
devices.  The  key  additional  feature  of  the  QHD  equations  is  the  incorporation  of  the  effects  of  density 
gradients  through  the  quantum  potential.  Insofar  as  a  classical  solution  does  not  exist  in  the  presence 
of  barriers,  Q  (i.e.,  one  of  die  major  quantum  correction  terms)  must  at  least  partially  caned  die  effects 
of  the  barriers  and  permit  carrier  transport  Indeed,  as  is  shown  in  Fig.  l,  Q/i  does  approximately  bal¬ 
ance  die  barrier  potential  which  is  a  result  consistent  with  Ref.  IS.  hi  this  work,  the  quantum 
corrected  HTM  was  used  to  demonstrate  that  die  hydrodynamic  approach  can  be  used  to  study 
ultrasmall  dectron  devices. 

In  another  simulation  study,  a  simplified  form  of  the  multi-valley  hydrodynamic  transport  model 
was  used  to  investigate  the  feasibility  of  utilizing  the  vdocity  overshoot  effect  in  reducing  electron 
transit  time  through  submicron  GaAs  electronic  device  structures.  In  this  initial  investigation  with  die 
new  multi-valley  model,  which  incorporates  a  separate  set  of  hydrodynamic  equations  for  each  conduc¬ 
tion  band  (T,  L,  and  X)  and  allows  for  electron  exchange  through  energy  dependent  relaxation  parame¬ 
ters,  the  nonparabobc  correction  terms  were  be  suppressed.  This  was  done  to  simplify  the  analysis  and 
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Figure  1:  Conduction  band  (£«)  and  quantum  potential  (Q)  at  =  8mF.  The  source 

and  collector  regions  are  excluded  from  the  plot. 
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to  provide  a  base  for  future  nonparabolic  studies.  This  new  multi-valley  model  was  applied  to  a  novel 
submicron  GaAs  ballistic  diode  structure  with  bandgap  engineered  electric  field  spikes  (and/or  doping 
spikes)  in  older  to  reduce  upper  valley  occupancy  and  enhance  device  transit  time  using  the  velocity 
overshoot  effect  The  general  results  of  the  study  (summarised  in  Fig.  2)  indicate  that,  while  velocity 
speed-up  is  minimal  in  the  neighborhood  of  the  active  (low-doped  or  high-electric  field)  region,  velo¬ 
city  overshoot  can  be  extended  into  the  collector  (low-electric  field)  region  with  significant  transit  time 
enhancement  Since  these  results  are  consistent  with  other  independent  investigations  using  more  accu¬ 
rate  (Monte  Carlo)  techniques,  this  work  has  proven  that  HTMs  are  a  viable  alternative  to  these  more 
computationally  intensive  methods. 

We  have  also  determined  the  impact  of  k-space  transfer  (Le.,  intervalley  transfer)  and  band  non- 
parabolicity  on  electron  were  determined.  In  this  work,  variations  of  the  multi-valley  (T,  L,  and  X) 
nonparabolic  hydrodynamic  transport  model  was  used  to  study  a  submicron  GaAs  ballistic  diode. 
Numerical  simulations  indicate  that  accurately  including  the  effects  of  nonparabolicity  in  the  streaming 
terms  and  k-space  transfer  in  tile  velocity  and  energy  equations  is  very  important  in  correctly  determin¬ 
ing  the  conductance  of  the  device.  The  existence  and  amount  of  negative  differential  conductance  was 
determined  to  be  strongly  influenced  by  both  of  these  physical  factors.  Furthermore,  tire  sensitivity  of 
device  conductance  to  changes  in  the  thermal  conductivity  is  diminished  significantly  when  nonparabol¬ 
icity  is  accurately  incorporated.  Results  for  the  total  current  density  versus  applied  bias,  for  different 
variations  in  the  model  are  given  in  Fig.  3.  These  simulations  have  demonstrated  the  pronounced 
effects  of  conduction  band  nonparabolicity  and  k-space  transfer  on  electron  transport  in  a  GaAs  ballis¬ 
tic  diode.  Therefore,  the  study  of  nonequilibrium  and  nonstationary  electron  physics  should  be  per¬ 
formed  with  this  type  of  advanced  macroscopic  transport  model 

More  recently,  work  has  been  performed  to  develop  an  extended  Scharfetter-Gummel  current- 
density  discretization  formula  sufficient  for  solving  an  advanced  hydrodynamic  electron  transport  model 
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Figure  2:  Outline  of  transport  study  and  transit  time  results 
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Figure  3:  Total  electron  current  density  Terras  applied  electric  potential  for  a  subxnicran 
ballistic  diode  (BD  sac  inset).  Results  are  for  variations  in  the  hydrodynamic  transport 
model  (HTM  ;  dashed  lines)  and  a  Monte  Carlo  study  (MC  ;  solid  line).  Four  different 
variations  (cases)  fat  the  HTM  were  consider.  The  results  for  these  different  cases  are 
designated  hy  dashed  lines  with  dash  lengths  that  increase  fat  the  following  order.  Case  A 
:  HTM  without  nanparabolidty  (in  the  streaming  terms)  or  ftU  fattervalley  effects  (in  the 
velocity  or  energy  equation).  Case  B :  HTM  with  fattervalley  effects  included  only.  Case  C : 
HTM  with  nonparabolicity  effects  included  only.  Case  D  :  HTM  with  nanparabolidty  and 
fattervalley  effects. 
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(HTM).  While  this  extended  formulation  property  ineoipontes  many  advanced  effects  associated  with 
hot-electron  transport,  its  primary  feature  is  the  accurate  treatment  of  the  convection  term  (or  displace¬ 
ment  energy  term)  which  is  present  in  both  conventional  (i.e.,  drifted  Maxwellian  model)  and  the 
advanced  nonparabolic  model  previously  discussed  in  this  report  In  general,  since  die  success  of  any 
numerical  simulation  (for  a  wide  variety  of  device  structures  and  operating  conditions)  is  critically 
dependent  on  convergence  rate  and  numerical  stability,  the  choice  of  an  efficient  discretization  scheme 
is  very  important  Currently,  most  robust  discretization  formulations  for  HTMs  utilize  extensions  of 
the  Scharfetter-Gummel  (SG)  scheme  which  was  found  to  be  stable  in  drift-diffusion  (DD)  model 
simulations.  However,  most  extensions  of  the  SG  formula  for  HTMs  either  ignore  convention  terms  or 
else  assume  a  priori  feat  they  are  defined  in  a  previous  step  of  an  iteration  procedure.  A  simple 
analytical  approach  for  directly  including  convection  effects  into  an  extended  SG  formulation  has  been 
established  and  several  versions  of  this  discretization  are  under  investigation.  Presently,  foe  stability 
and  efficiency  of  this  new  discretization  formulation  is  still  under  investigation. 

Finally,  work  has  been  performed  to  interrogate  foe  influence  of  different  device  structures  and 
boundary  conditions  on  the  resultant  heat  flow  vector  (a  transport  parameter  in  the  HIM)  using  the 
Monte  Carlo  Approach.  Thus  far,  we  have  observed  that  die  sign  of  the  heat  flow  vector  is  a  function 
of  position  in  ballistic  diode  structures  and  that  boundary  conditions  do  influence  the  effective  conver¬ 
gence  rate  (i.e.,  simulation  time  to  achieve  correct  symmetry  of  file  electron  distribution)  of  file  Monte 
Carlo  method. 

2.2.2  Modeling  and  characterization  of  novel  high  speed  devices 

During  die  past  contract  period,  we  have  studied  electron  transport  properties  and  electrical 
characteristics  for  a  number  of  high  speed  device  structures,  including  delta-doped  GaAs  MESFETs, 
AlInAs/GalnAs/InP  HEMTs  with  different  additional  channel  doping  characteristics,  novel  real-space 
transfer  logic  transistor  structures,  and  devices  with  variable  doping  profiles.  The  motivations  of  our 
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research  on  these  structures  are  to  provide  in-depth  understanding  of  electron  transport  in  devices,  to 
investigate  the  possibility  of  realizing  novel  high  speed  functional  device  structures,  and  to  provide 
general  design  guidelines  for  optimal  device  operation. 

As  a  promising  high  speed  device,  the  delta-doped  OaAs  MESFET  features  high  current  drive 
capability,  high  transconductanoc.  improved  threshold  voltage  control,  and  improved  breakdown 
characteristics.  The  study  of  properties  of  delta-doping  in  OaAs  predicts  that  device  performance  of 
delta-doped  OaAs  MESFETs  is  comparable  or  even  better  than  those  of  AlGaAs/GaAs  HEMTs  [19]. 
However,  experimental  results  show  a  wide  diversity  and  fell  short  of  theoretical  predictions.  The 
design  and  optimization  of  device  structure  (especially  delta-doping  profile)  play  an  important  role  in 
guiding  die  improvement  of  device  performance.  Another  high  speed  device,  the  AlInAs/GalnAs/InP- 
based  HEMT,  has  been  the  focus  of  extensive  experimental  study  due  to  superior  electron  transport 
properties  in  this  lattice-matched  material  system.  While  promising  dc  and  microwave  performances 
have  also  been  achieved.  Ugh  output  conductance  has  been  repeatedly  observed  for  devices  with 
channel-length  in  die  deep-submicron  (less  than  0.3  pm)  regime  [20,21].  Physics  insight  leading  to  the 
Ugh  output  conductance  and  potential  solutions  merit  careful  investigation  if  potential  advantages  of 
AlInAs/GalnAs  HEMTs  are  to  be  fully  exploited.  The  third  device  structure  which  we  have  investi¬ 
gated  this  year  belongs  to  a  new  class  of  devices  which  employs  die  concept  of  electron  real-space 
transfer  [22]  in  heterojunction  structures.  The  charge  injection  transistor  (CHINT)  [23]  and  the 
recently-developed  CHINT  logic  element  [24]  are  representative  devices.  Theoretical  predictions  and 
experimental  evidence  from  CHINT  operation  indicate  the  possibility  for  developing  novel  and 
improved  multi-terminal  real-space  transfer  devices  which  can  realize  multiple  logic  functions.  Our 
last  work  in  the  area  of  device  modeling  deals  with  the  effects  of  electron  transport  in  devices  with 
variable  doping  distributions.  Such  device  structures  can  be  realized  by  utilizing  advanced  fabrication 
techniques,  such  as  focused  ion  beans  or  epitaxy  (MBE,  MOCVD,  ALE)  methods.  The  potential 
advantages  of  variable  doping  in  device  channel  or  active  transport  region  can  best  be  realized  when 


detailed  electron  transport  properties  in  practical  device  structures  are  understood  and  doping  scheme 
optimization  is  achieved.  In  die  following,  we  outline  the  main  results  of  our  theoretical  study.  More 
detailed  descriptions  for  each  subject  can  be  found  in  the  related  papers  and  presentations. 

A.  Delta-doped  GaAs  MESFET 

Figure  4  gives  a  schematic  illustration  of  the  device  structure  which  we  have  simulated,  along 
with  the  delta-doping  profile.  In  this  work,  we  use  a  twodimensional,  drift-diffusion  model  to  exam¬ 
ine  die  influence  of  key  design  parameters  on  device  operation.  These  parameters  are  (Fig.  4)  peak 
doping  density,  half  width  at  half  maximum  of  the  doping  profile,  and  distance  horn  the  location  of 
peak-doping  to  die  device  surface,  top  layer  background  doping,  and  lateral  feature  sizes.  The  general 
design  considerations  and  trade-ofis  are  studied.  In  particular,  our  simulation  results  reveal  than  1)  A 
steeper  doping  profile  and  a  short  gate-to-channel  distance  are  preferred  for  improvement  in  device  per¬ 
formance.  2)  The  peak  delta-doping  concentration  and  half-width  at  half  maximum  of  delta-doping 
profile  are  found  to  have  a  more  profound  effect  than  does  delta-doping  depth  (dp)  in  enhanced  overall 
device  performance.  And  3)  careful  trade-off  should  be  made  in  the  design  of  top-layer  background 
doping  in  terms  of  current  drive,  transconductance,  and  threshold  voltage  sensitivity  control  Our 
results  show  that  a  relatively  high  top  layer  background  doping  is  preferred  to  increase  drive  current 
and  reduce  relative  threshold  voltage  shift  (due  to  processing  parameter  fluctuations)  at  a  cost  of  a 
modest  reduction  in  transconductance.  Simulation  results  are  in  dose  agreement  with  measurements 
horn  experimental  devices  fabricated  at  Texas  Instruments  using  atomic  layer  epitaxy  material  prepared 
in  NCSU. 

B.  AUnAslGaJnAs  HEMT 

Figure  5  shows  the  simulated  quarter-micron-gate  AlInAsA3alnAs  HEMT  structure  which 
employs  different  additional  charnel  doping  (h\  i,  and  p)  configurations.  The  motivation  is  to  study 


Figure  4:  (a)  mnstretion  of  a  delta-doped  GaAs  MESFET.  The  mesh  distribution  is  illus¬ 

trated  in  the  imeL  b)  Delta-doping  profile.  The  solid  line  shows  the  measurement  data  and  the  dashed 
line  shows  the  distribution  used  in  the  simulations. 
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Figure  5:  drswing  of  the  layer  structures  for  AlInAs/GalnAs  HEMTs  used  in  this  study. 
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die  effect  of  additional  channel  doping  on  device  transconductance,  output  conductance,  and  high  fre¬ 
quency  performance.  The  AlInAs/GalnAs/InP  material  system  appears  to  be  the  most  promising  candi¬ 
date  for  high  frequency  applications.  However,  for  this  material  system,  high  output  conductance  has 
been  an  major  problem  that  degrades  high  frequency  performance  when  the  HEMT  device  feature  size 
is  reduced  to  the  ultrasubmicron  (less  than  0.3  pm)  regime.  This  has  been  reported  by  several  research 
groups.  Attempts  have  bear  made  experimentally  to  reduce  output  conductance,  however,  physical 
insight  remains  an  unknown.  In  this  work,  we  investigated  the  possible  mechanisms  which  lead  to 
enhanced  or  degraded  device  performance  in  different  channel -doped  structures  using  two-dimensional, 
self-consistent  ensemble  Monte  Carlo  simulations.  Our  results  indicate  that  conduction  in  the  AlInAs 
buffer  layer  and  reduced  electric  field  shielding  in  n-channel  doped  HEMT  are  the  possible  reasons  for 
increased  output  conductance.  By  use  of  a  properly  designed  p-doped  channel  layer,  we  observe 
significantly  reduced  output  conductance  and  increased  transconductance  compared  with  those  of  n- 
channel  doped  HEMTs.  Detailed  analysis  indicates  that  improved  electron  confinement  in  the  device 
channel  and  increased  field  shieding  for  HEMT  with  p-channel  doping  improve  overall  device  perfor¬ 
mance.  Our  simulations  are  in  reasonable  agreement  with  measurement  results  from  experimentally 
fabricated  devices  [with  different  additional  channel  doping  (UCSB  and  Hughes  Research  Lab)]. 

C.  Real-space  transfer  logic  transistor 

The  basic  teal-space  transfer  logic  transistor  (RSTLT)  is  a  four-terminal  device,  which  has  two 
separate  collectors  (Cl  and  C2)  as  output  terminals  when  connected  through  load  resistors  (Rci  and 
Res)  to  the  collector  supply  voltage  (Vc);  a  drain  input  terminal  (D);  and  a  common  source  terminal 
(S)  which  is  at  ground  potential  (Fig.  6).  The  operational  principle  of  the  proposed  RSTLT  is  based 
on  the  localized  nature  of  real-space  transfer  of  hot  electrons  in  such  structures.  Channel  electrons 
gain  energy  from  the  lateral  electric  field  and  real-space  transfer  occurs  when  foe  energy  of  electrons  is 
comparable  to  or  higher  than  the  heterointerface  barrier  height  The  transferred  electrons  will  be  col- 
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lected  by  either  Cl  or  C2  through  vertical  electric  field.  For  a  given  heterointerface  barrier  height,  the 
spatial  location  of  electron  RST  is  determined  primarily  by  the  bias  conditions.  By  properly  designing 
the  device  configuration  and  selecting  die  operating  voltages,  ultra-fast  current  exchange  between  two 
collectors  can  be  controlled  by  the  applied  drain  voltage.  Compared  with  CHINT  structure,  the  pro¬ 
posed  RSTLT  employs:  1)  a  recessed  Cl  collector  to  enhance  die  vertical  electric  field  intensity  under 
Cl  for  matched  output  collector  currents,  2)  unintentionally  doped  channel  and  barrier  layers  for 
reduced  collector-source  and  collector-drain  leakage  currents,  and  3)  a  collector-up  device  configuration 
to  reduce  parasitic  capacitance  and  to  facilitate  proper  arrangement  of  the  device  terminals. 

In  this  project,  we  employ  a  self-consistent  ensemble  Monte  Carlo  model  to  perform  die  device 
feasibility  study.  Figures  7  and  8  show  die  current-voltage  characteristics  and  transient  current 
response  as  a  function  of  time,  respectively.  A  characteristic  delay  time  of  about  3.0  psec  can  be 
achieved,  which  would  make  the  RSTLT  one  of  the  fastest  switching  devices  reported  to  dale.  This 
characteristic  current  switching  (between  two  collector  terminals)  can  be  utilized  to  realize  a  multi-logic 
element  in  a  single  device  structure.  We  have  performed  case  studies  to  demonstrate  that  by  properly 
selecting  device  dimensions  and  the  resistor  load  elements,  collector  terminals  Cl  and  C2  can  function 
as  logic  NOT  and  EQUIVALENT  of  input  terminal  D.  Based  on  these  results,  we  propose  a  RSTLT 
logic  element  (Fig.  9),  which  features  NORIAND  functions  (Cl  functions  as  NOR  of  D1  and  D2;  C2 
functions  as  AND  of  D1  and  D2,  respectively)  and  a  comparable  characteristic  delay  time  as  that  of  die 
basic  RSTLT. 

D.  Effects  of  doping  variation  on  electron  transport 

High  performance  electronic  systems  require  devices  with  high  current  drive  capability,  high  tran¬ 
sconductance,  high  breakdown  voltage,  and  low  parasitic  resistance.  Currently,  most  semiconductor 
devices  employ  doping  techniques  that  control  doping  profiles  in  die  vertical  dimension  (e.g.,  uniform 
doping,  modulation  doping,  and  delta-doping),  that  is,  perpendicular  to  the  direction  of  carrier 


Current  (mA/mm) 


Figure  8:  Current  response  for  drain  voltage  switch-off  for  an  AlojGaQjAs/GaAs  RSTLT.  The  col¬ 
lector  voltages  are  fixed  at  Vci«Vqs4.0  V.  The  device  input  (Vjj)  is  switched  from  4  V  to  1  V  at 
t*=1.0  psec. 
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transport  The  rapid  development  of  new  technologies  makes  possible  changes  in  doping  profile  in 
both  lateral  and  vertical  (parallel  and  normal  to  carrier  transport  respectively).  This  added  degree  of 
freedom  for  "doping  engineering"  can  possibly  result  in  novel  and  improved  device  structures.  The 
motivations  of  this  project  are  to  study  die  effects  of  unique  variable  doping  profiles  on  electron  tran¬ 
sport  properties  and  to  assess  potential  device  applications. 

We  apply  self-consistent  ensemble  Monte  Carlo  model  to  study  a  GaAs  n+-n~n+  test  device. 
Doping  schemes  are  ramp-doping  (ramp-up  and  ramp-down)  and  spike -doping  (single-spike  and  multi - 
spike)  in  the  n  region.  In  order  to  have  meaningful  comparisons  among  these  structures,  we  apply 
conditions  of  identical  minimum  doping  density  and  the  same  total  integrated  dopant  in  the  active  n 
region.  Figure  10  shows  the  current-voltage  characteristics  of  ramp-up  (from  the  left  n+-n  junction), 
ramp-down,  and  two  constant  doped  structures.  Figure  1 1  compares  the  current-voltage  characteristics 
for  one-spike  and  three-spike  doped  structures.  It  is  found  for  conditions  of  identical  minimum  doping 
density  and  total  integrated  dopant  in  the  n  region  that:  1)  the  most  significant  improvement  can  be 
achieved  by  employing  a  structure  with  a  ramp-down  doping  scheme.  2)  Introducing  n+  spike(s)  in  the 
active  electron  transport  (n)  region  brings  improved  device  performance  compared  to  devices  without 
n+  spikes.  For  the  one-spike  doped  structure,  electron  transport  is  independent  on  the  location  of  the 
n+  spike.  However,  less  significant  improvement  is  observed  for  this  structure  due  to  the  presence  of 
retarding  fields  over  a  large  portion  of  the  n  region.  3)  Further  improvement  can  be  achieved  by  split¬ 
ting  one  n+  spike  into  several  spikes  in  the  n  region.  The  redistribution  of  electron  density  and  the 
reduced  effect  of  retarding  electric  fields  bring  an  increased  average  electron  velocity  and  a  resulting 
higher  current  drive  capability.  These  results  indicate  the  possibility  of  enhanced  device  performances 
for  FET  devices  which  employ  proper  variable  doping  schemes.  Potential  advantages  of  variable  dop¬ 
ing  in  normally-on  FET  device  applications  include  enhanced  current  drive,  reduced  source  resistance, 
and  improved  breakdown  characteristics. 


Current  (106  A/cm2) 


Voltage  (V) 


Figure  10:  Current-voltage  characteristics  for  ramp-doped  n+-n— n+  structures.  The  doping  density 
varies  linearly  between  5xl016  cm"3  and  5xl017  cm*3.  Simulation  results  for  n+-n-n+  structures  with 
constant  n  doping  densities  of  5xl016  cm*3  and  1.6xl017  cm*3  are  included  for  comparison. 
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Figure  11:  Current-voltage  characteristics  for  one-spike,  three-spike,  and  constant  (SxlO16  cm"3) 
doped  n+-n-n+  structures.  Hie  doping  densities  of  the  n+  spike(s)  and  the  background  are  9.16x10" 
cm"3  and  SxlO16  cm'3,  respectively.  The  one-spike  structure  has  die  n+-spike  placed  in  die  center  of 
the  n  region. 
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2.2.3  Modeling  of  hot  elections  in  silicon  devices 

For  the  past  fifteen  years,  the  silicon  MOSFET  industry  has  been  dealing  increasingly  with  prob¬ 
lems  related  to  hot  electron  injection  intc  .*  oxide  with  production  mode  devices.  For  even  longer, 
laboratories  have  been  measuring  and  modeling  hot  electron  injection  into  die  oxides  with  various  lev¬ 
els  of  success  [25].  Frequently,  the  model  of  choice  has  been  the  "lucky  electron”  model  coupled  with 
a  drift-diffusion  simulator.  While  the  "lucky  electron"  model  has  enjoyed  a  reasonable  amount  of  suc¬ 
cess,  we  have  recently  entered  the  regime  [26]  where  drift-diffusion  simulation  is  no  longer  accurate. 
With  this  reality,  researchers  have  increasingly  been  forced  to  switch  to  hydrodynamic  and  Monte 
Carlo  simulators  to  better  predict  the  behavior  of  devices.  Using  these  models  (especially  Monte 
Carlo),  significantly  more  information  about  the  electron  energy  distribution  is  available,  giving  us  fire 
opportunity  to  improve  upon  the  "lucky  electron"  model.  Full  theoretical  treatment  of  die  entire  hot 
electron  injection  condition  poses  many  serious  technical  difficulties  [27]  but  some  first  pass  efforts 
[28]  indicate  we  may  have  a  reasonable  hope  for  success.  We  can  then  couple  these  injection  charac¬ 
teristics  with  the  oxide  trapping  characteristics  [29,30]  to  give  us  a  second  generation  model  to  sup¬ 
plant  the  "lucky  electron"  model. 

Within  the  past  year,  we  have  developed  a  state-of-the-art  Monte  Carlo  simulator  of  electrons  for 
silicon  devices.  This  simulator  uses  a  realistic  band  structure  (derived  from  pseudopotential  calcula¬ 
tions)  for  die  best  accuracy  in  simulating  band  effects  on  electron  simulation.  It  also  uses  an  ensemble 
method  of  modeling  which  allows  the  simulation  to  account  for  both  spatial  and  electron-electron  inter- 
collisional  effects  in  a  self-consistent  manner.  In  addition,  our  Monte  Carlo  model  includes  all  of 
major  scattering  mechanisms  important  for  silicon  simulation.  These  mechanisms  are  acoustic  phonon, 
intervalley  phonon,  impact  ionization,  ionized  impurity  and  electron-electron  scatterings.  This  rather 
accurate  implementation  of  physical  details  results  in  a  reasonable  field  vs.  velocity  and  scattering  rates 
which  compare  favorably  with  die  full  numerical  calculations  of  Fischetti  [31]. 
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In  order  to  accurately  model  complex  device  structures,  we  have  also  developed  a  non-uniform 
mesh  Poisson  solver  which  is  re-evaluated  at  frequent  time  intervals  to  reflea  changes  in  the  potential 
distribution  due  to  the  movement  of  the  electrons  in  die  simulation.  This  movement  is  mirrored  in  the 
calculation  of  background  hole  concentrations.  The  carrier  concentrations  remain  self-consistent  by 
updating  the  hole  concentrations  using  the  quasi-Fermi  level  derived  from  die  Poisson  solution.  At  die 
same  time,  an  efficient  method  has  been  implemented  which  allows  to  transfer  PISCES-modeled  and 
developed  structures  into  one  of  our  Monte  Carlo  runs.  This  approach  will  improve  overall  computing 
time  and  efficiency  significantly  by  providing  a  smart  initial  guess  and,  thus,  will  permit  us  to  concen¬ 
trate  to  the  regimes  of  sensitivity  in  a  greater  detail. 

Calculations  of  actual  MOSFET  structures  is  now  finishing  the  debugging  stage  and  the  prelim¬ 
inary  data  begin  to  show  the  importance  and  value  of  our  model  in  die  study  of  hot  electron  effects 
and  the  related  aging  issues  in  these  devices.  Figure  12  demonstrates  that  elections  experience  a 
significant  velocity  overshoot  region  even  at  reasonably  long  channels  (l  pm)  when  stressed  at  high 
voltages  (V*=6  V,  Vf=3  V,  V,=Vb=0  V).  The  electron  energy  distribution  tail  in  die  drain  pinch-off 
region  is  also  seen  to  have  a  significant  population  beyond  die  3.1  eV  necessary  to  inject  electrons  into 
the  oxide  as  shown  in  Fig.  13.  Such  details  are  the  essential  information  to  understand  nonlinear  elec¬ 
tron  transport  in  submicron  devices.  We  expea  to  obtain  more  data  in  our  effort  to  develop  a  realistic 
model  for  oxide  trapping  and  die  generation  of  interface  states  at  the  Si/Si  C>2  interface  within  die  next 
few  weeks. 

23  Publications  and  Presentations 

During  the  last  year  we  have  made  six  oral  presentations  at  national  and  international  confer¬ 
ences.  Also,  twelve  written  manuscripts  have  been  published  in  the  refereed  literature.  In  addition, 
another  manuscript  has  been  accepted  for  publication  and  six  manuscripts  have  been  submitted  or  are 
currently  in  preparation.  1116  following  paragraphs  summarize  die  presentations  and  publications  made 
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Figure  12:  Calculated 
electron  velocity  along 
the  channel  in  a  1-ym 
gate-length  MOSFET  at 
higli  voltages  (Vf=6  V, 
V,=3  V,  V^V*=0  V). 


Figure  13:  Calculated 
election  energy  distri¬ 
bution  along  the  chan¬ 
nel  for  the  same  device 
and  bias  conditions  as 
in  Fig.  12. 
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under  this  program  during  the  last  year,  and  include  material  not  previously  reported  to  ONR  in  the 
last  annual  report  as  well  as  material  which  was  accepted  for  publication  in  1991  and  has  since  been 
published  in  1992. 
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We  present  results  from  a  theoretical  study  of  the  influence  of  doping  profile  variations  on 
the  performance  of  delta-doped  AIGaAs/GaAs  high  electron  mobility  transistors 
(HEMTs).  An  ensemble  Monte  Carlo  simulation  coupled  with  a  self-consistent  solution  of 
the  two-dimensional  Poisson  equation  is  used  to  investigate  HEMTs  which  employ 
both  single  and  double  delta-doped  profiles  with  varying  doping  configurations.  The  calculated 
results  reveal  that  single  delta-doped  HEMTs  designed  with  identical  threshold  voltages 
exhibit  improved  device  behavior  when  thinner  delta-doped  layers  with  more  heavily  doped 
concentrations  an  utilized.  For  double  delta-doped  HEMTs  with  an  identical  total 
doping  in  the  AlGaAs  layer,  improved  threshold  voltage  control  is  obtained  as  the  spacing 
between  two  delta-doped  layers  increases.  However,  this  increase  in  spacing  also 
causes  a  degradation  in  transconductance,  cut-off  frequency,  and  switching  time.  As  gate  bias 
increases,  the  dependence  of  device  performance  (or  degradation)  on  the  spacing 
between  doping  planes  becomes  less  pronounced  due  to  the  upward  shift  in  threshold  or 
“onset”  of  parallel  conduction  in  the  AlGaAs  layer. 


I.  INTRODUCTION 

The  high  electron  mobility  transistor  (HEMT)  was 
first  developed  in  19801  and,  today,  the  AIGaAs/GaAs 
HEMT  remains  as  one  of  the  most  important  heterojunc¬ 
tion  structures  due  to  its  excellent  device  performance,  the 
maturity  of  the  materials  growth  technology,  and  the  ex¬ 
cellent  lattice  match  for  the  AIGaAs/GaAs  material  sys¬ 
tem.  Further  device  improvements  are  needed  in  order  to 
fully  exploit  the  potential  of  AIGaAs/GaAs  HEMT  struc¬ 
tures.  Unfortunately,  improvements  in  “conventional”  uni¬ 
formly  doped  HEMTs  are  limited  by  the  occurrence  of 
persistent  photoconductivity,1  threshold  voltage  shift,1  and 
collapse  of  current-voltage  (I-V)  characteristics4  due 
largely  to  effects  caused  by  DX  centers  and  surface  states. 
The  use  of  low  A1  mole  fraction  can  reduce  these  effects; 
however,  this  also  reduces  the  conduction-band  edge  dis¬ 
continuity,  resulting  in  decreased  channel  electron  density 
and  degraded  device  performance.  This  inherent  limitation 
can  be  circumvented  by  device  structure  modifications, 
such  as  those  employed  in  delta-doped  (also  referred  to  as 
pulse-doped  or  planar-doped)  AIGaAs/GaAs  HEMT 
structures.  Recently,  delta-doped  HEMTs  have  attracted 
much  attention.  Experiments  show  that  these  devices 
achieve  improved  threshold  control,  high  channel  electron 
density,  reduced  trapping  effects,  and  high  breakdown 
characteristics.*-10 

Detailed  investigations  of  delta-doped  HEMTs  are 
necessary  to  reveal  and  understand  the  underlying  device 
physics.  While  extensive  study  of  conventional  uniformly 
doped  HEMTs  has  been  reported, IM>  the  theoretical  in¬ 
vestigation  of  delta-doped  HEMTs  is  rare.*'’’1'  The  opera¬ 
tion  of  submicron  delta-doped  HEMTs  involves  nonsta¬ 
tionary  electron  transport  and  real-space  transfer 
effects. ,MI  Also,  the  degeneracy  effect,  DX  centers  and 


surface  states  have  significant  influence  on  practical  device 
operation. 19-14  Therefore,  the  development  of  more  accu¬ 
rate  and  complete  device  models  must  include  these  effects. 
Furthermore,  it  is  important  to  explore  new  approaches 
which  can  lead  to  improvements  in  AIGaAs/GaAs 
HEMTs.  Specifically,  reducing  series  resistance  and  in¬ 
creasing  current  drive  capability  of  the  single  delta-doped 
HEMT  structure  is  of  importance  for  microwave  and 
power  applications.  As  an  alternative,  a  structure  with  a 
second  delta-doped  layer  introduced  in  the  AlGaAs  layer 
(i.e.,  a  double  delta-doped  HEMT)  could  be  used.  This 
structure,  while  unavailable  experimentally,  merits  study 
since  the  introduction  of  a  second  delta-doped  layer  could 
possibly  provide  more  drive  current  and  reduce  source-gate 
and  gate-drain  series  resistance  without  reducing  the  ad¬ 
vantages  provided  by  single  delta  doping.  A  double  delta- 
doped  structure  can  be  realized  by  molecular-beam  epitaxy 
or  metalorganic  chemical  vapor  deposition  technologies. 
However,  a  theoretical  study  can  provide  guidance  for  the 
experimentalists  to  assess  the  performance  of  such  a  device 
and  its  dependence  on  doping  parameters  and  bias  condi¬ 
tions. 

With  these  motivations,  we  employ  physical  device 
models  to  examine  the  effects  of  doping  profile  variations 
on  the  performance  of  two  sets  of  HEMT  revices.  The  first 
set  includes  single  delta-doped  HEMTs  with  varying  dop¬ 
ing  configurations  and  constant  threshold  voltage.  Aim, 
for  the  first  time,  an  attempt  is  made  to  analyze  the  per¬ 
formance  of  proposed  double  delta-doped  HEMTs.  This 
second  set  of  devices  is  studied  by  varying  the  spacing 
between  two  delta-doped  layers  with  identical  widths  and 
total  doping  concentrations.  Device  performance  is  exam¬ 
ined  in  terms  of  threshold  voltage  control,  transcondoc- 
tance,  cut-off  frequency,  and  characteristic  switching 
times.  [Since  the  main  purpose  of  this  work  is  to  study  the 
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A  pstb-snisyru!  Metric  Ceric  method  Car  calculation  of  real  self-energies  of  free 
and  conined  carriers  at  non  sere  temperatures  is  presented.  An  efficient  Burner- 
teal  method  fot  lh<  polii  optical  phoaoa  wftunct  ^wflww?  hat  been 

developed  and  incorporated  into  a  previously  developed  P1MC  for  calculating  equi¬ 
librium  properries  of  uncoupled  carriers.  As  illustrative  examples,  calculations  of 
carrier  teal  sclLeaergies  are  performed  for  bulk  0str  as  a  function  of  temperature, 
and  far  AIGaAs-GcAa-likc  square  quantum  vires  at  77K  and  300K  as  a  function  of 
wire  width.  The  self-energies  obtained  here,  in  part,  represent  new  results  far  these 
systems.  Where  previous  results  are  available,  the  agreement  is  good. 


1.  Introduction 

The  Aysmss  peihwiqis/  (FPl)  formulation  of  quan¬ 
tum  mechanics  allows  formal  inclusion  of  carrier  phonon 
coupling  to  all  orders  in  the  coupling  potential  |1,2].  For 
this  reason  the  FPI  formalism  long  has  been  used  to  ana¬ 
lytically  calculate  carrier  real  self-energies  and  effective 
masses  in  bulk  semiconductors  (8-4).  However,  exten¬ 
sion  of  the  FPI  formalism  to  analytic  study  of  ultrasmaO 
devices  has  been  hindered  by  approximations  required  to 
model  confined  geometries  |7|.  Previously,  though  im¬ 
proved  results  were  not  obtained,  it  was  demonstrated 
that  carrier  self  energies  in  bulk  semiconductors  at  9K 
could  he  calculated  using  numerical  Paik-lnltfrtl  Metric 
CsrU  (PIMC)  methods  as  well  (8-19).  However,  to  the 
best  of  the  authors’  knowledge,  these  calculations  were 
not  extended  to  aoaseto  temperatures  or  confined  do- 
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vice  geometries.  Here,  a  PIMC  method  far  calculating 
free  and  confined  carrier  self-energies  at  non-sero  tem¬ 
peratures  is  presented. 

2.  Feynman  Path-Integral  Formalism 

la  the  FPI  formalism,  for  a  coupled  carrier-phonon  sys¬ 
tem,  the  equilibrium  density  matrix  far  the  carrier  av¬ 
eraged  ever  the  phonon  coordinates  and  subject  to  a 
MaxweU-Boltsmann  distribution  in  energies  at  temper¬ 
ature  T  is  given  by 

MfAT)  -  j^*p{- J +  FROIW*} 

fl) 

Here,  the  *j£...Px(t)a  signifies  an  integral  ever  all 
paths  from  fi  to  Fthat  ate  continuous  functions  of  time, 
and  is  the  imf fames  /unctions!  that  accounts  far 
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This  paper  presents  a  hydrodynamic  model  suitable  for  studying  hot-electron  transport  in  semicon¬ 
ducting  materials  with  non  parabolic  conduction  bands.  The  model  presented  is  based  upon  a  unique 
derivation  of  the  moments  of  the  Boltzmann-transport  equation  for  the  streaming  (collision-independent) 
terms.  This  derivation  implements  an  efficient  and  compact  mathematical  formalism  appropriate  for 
electrons  under  the  influence  of  high  electric  fields  and  nonstationary  conditions.  The  theoretical  inves¬ 
tigation  also  introduces  a  distributional  form  with  nonparabolic  properties  to  precisely  define  the  result¬ 
ing  non  parabolic  streaming  parameters.  The  final  set  of  model  equations  is  exhibited  in  a  fashion  to 
clearly  show  the  correction  factors  to  the  more  familiar  hydrodynamic  model  applicable  for  the 
constant-eflective-mass  case.  In  general,  the  hydrodynamic  (or  conservation)  model  contains  pure  trans¬ 
port  terms  that  are  treated  as  being  independent  of  the  specific  dissipation  mechanisms  and  collision 
terms  to  directly  account  for  the  influence  of  scattering.  Since  the  collision  terms  are  almost  always 
treated  phenomenologically  using  a  relaxation-time  approximation,  our  formulation  of  the  streaming 
terms  should  significantly  improve  the  overall  accuracy  of  the  approach.  In  addition,  this  paper  presents 
the  results  of  an  extensive  investigation  of  the  assumed  ansatz  distribution  and  resulting  nonparaboiic- 
model  parameters  using  an  elaborate  Monte  Carlo  model.  The  Monte  Carlo  technique  was  used  to  gen¬ 
erate  comparison  electron  distributions  and  exact  values  for  the  non  parabolic  transport  parameters  for 
stationary  and  nonstationary  electronic  structures.  In  all  cases,  excellent  agreement  was  found  between 
the  Monte  Carlo-calculated  parameters  and  the  derived  non  parabolic-model  terms.  The  Monte  Carlo 
calculations  also  revealed  that  the  ansatz  distribution  used  in  the  derivation  represented  a  significant  im¬ 
provement  over  the  more  familar  displaced  Maxwellian.  Therefore,  this  model  should  prove  very  valu¬ 
able  for  studying  electronic -device  structures  operating  under  high-bias  conditions. 


I.  INTRODUCTION 

During  the  past  two  decades  there  has  been  an  in¬ 
creased  use  of  hydrodynamic  conservation  models  to  in¬ 
vestigate  nonstationary  and  nonequilibrium  electron  dy¬ 
namics  in  submicrometer  semiconductor  devices.1-14 
These  models,  referred  to  by  some1514  as  the  hydro¬ 
dynamic  equations  due  to  their  similarity  to  the  Enler 
equations  of  fluid  dynamics17,1*  used  in  classical  hydro¬ 
dynamics  studies,  are  based  upon  higher  moments  of  the 
Boltzmann-transport  equation  (BTE).  The  popularity  of 
the  hydrodynamic  electron  transport  theory  is  due  to  the 
physical  and  practical  attributes  of  the  approach.  The 
hydrodynamic  models  have  the  capability  to  include  non¬ 
stationary  and  hot-electron  effects'*  and  hence  are  superi¬ 
or  to  simple  drift-diffusion  (extreme  thermal  equilibrium 
approximation30  and  local-eiectric-field-dependent  case) 
models.  While  Monte  Carlo  methods  that  solve  the  full 
BTE  can  easily  incorporate  complicated  band  structures 
and  detailed  scattering  rates,  hydrodynamic  models  re¬ 
quire  much  less  computation  time  to  generate  solutions 
and  possess  macroscopic  terms  which  offer  important 
physical  insight  Thus  the  hydrodynamic  approach  offers 
much  flexibility  for  future  theoretical  electron  device 
studies. 

Stratton10  introduced  the  general  conservation  or 
momentum-energy  balance  approach  to  investigate  hot- 
electron  transport  in  semiconductors.  The  analysts  per¬ 


formed  by  Stratton  utilized  a  spherical  harmonic  expan¬ 
sion  with  the  relaxation-time  approximation  to  define  a 
nonequilibrium  distribution  which  was  only  slightly  an¬ 
isotropic.  This  analysis  replaced  the  nonequilibrium  dis¬ 
tribution  by  a  Maxwellian  distribution  function  in  some 
terms  of  the  BTE.  Thus  some  of  the  streaming  terms  in 
the  final  transport  equations  were  suppressed.  Blo- 
tefcjaer11  extended  the  theory,  using  a  more  general 
analysis  which  retained  all  terms  of  the  moment  equa¬ 
tions,  and  derived  relations  applicable  to  a  two-valley 
semiconductor.  These  previous  investigations  were  appl¬ 
icable  to  semiclassical  transport  and  assumed  homogene¬ 
ous  materials  with  parabolic  conduction  bands. 

Since  the  early  work  of  Stratton  and  of  Blotekjaer,  the 
hydrodynamic  approach  has  been  utilized  extensively  to 
study  various  transport  phenomena  in  many  different 
Si,*'*’22  Oa As, 1,13,14  and  heterostnicture11, 35,14  electron 
devices.  The  particular  material,  device  configuration, 
and  bias  conditions  under  consideration  have  led  to  a 
variety  of  acceptable  simplifications  and  to  a  hierarchy  of 
approximate  electron-transport  models  (for  example,  see 
Sandbom,  Rao,  and  Blakely17  and  references  therein). 
Abo,  much  effort  has  been  directed  toward  developing 
numerical  techniques  to  solve  the  hydrodynamic  trans¬ 
port  model  self-consistently  with  Poisson’s  equation  for 
many  semiconductor  structures.  This  has  included 
methods  for  solving  the  problem  with  and  without  time 
dependency  and  in  multiple  space  dimensions3*-30  as  well 
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Abstract 

A  new  multi-valley  hydrodynamic  transport  model  has  been  used  to  investigate  electron  transit  time 
through  eubmicroa  QaAs  electron  device  structures.  The  study  was  performed  to  determine  the  feasibility 
of  utilising  the  velocity  overshoot  effect  in  reducing  the  overall  transit  time.  This  new  model,  which 
incorporates  a  separate  set  of  hydrodynamic  equations  for  each  conduction  band  valley  (r,  L  and  A') 
and  allows  for  electron  exchange  through  energy  dependent  relaxation  parameters,  is  applied  to  a  novel 
structure  which  reduces  upper  valley  occupancy  and  enhances  device  transit  time.  The  results  of  this 
study  indicate  that  while  velocity  speed-up  is  minimal  in  the  neighborhood  of  the  active  (low-doped) 
region  that  the  velocity  overshoot  can  be  extended  into  the  collector  region  (high-doped)  region  with 
significant  transit  time  enhancement. 


Introduction 

In  recent  years,  a  great  interest  has  developed  in  utilising  velocity  overshoot  in  both  GaAs  homojuaction 
and  heterojunctioa  BJT’s  to  achieve  lower  electron  transit  times  and  hence  higher  operating  frequencies!  1, 
2,  3].  Investigations  of  these  devices  have  shown  the  emitter-collector  transit  time  (ne)  and  the  collector 
transit  time  (ra)  to  be  dominant  (actors  in  determining  the  total  transit  time(4].  This  fact  has  directed 
much  energy  in  the  simulation  and  modeling  communities  toward  designing  and  studying  electronic 
structures  which  would  reduce  these  figures  of  merit  (r„  and  r„).  In  fact,  there  have  been  many  theoretical 
results  (8, 6,  7]  which  indicate  that  velocity  overshoot  can  be  exploited,  by  select  designs  and  proper  biases, 
to  achieve  very  efficient  high  frequency  devices.  However,  questions  have  been  raised  as  to  whether 
conventional  device  structures  can  exhibit  such  enhanced  effects.  Also,  simulation  evidence  has  been 
reported  to  indicate  that  velocity  overshoot  may  be  over  estimated  in  such  structurcs(8]. 

In  this  paper,  a  new  multi- valley  hydrodynamic  transport  model  is  used  to  investigate  the  feasibility  of 
utilising  the  velocity  overshoot  effect  in  reducing  electron  transit  time  through  submicron  GaAs  electronic 
device  structures.  In  previous  work  (B,  10],  a  novel  approach  has  been-presented  to  incorporate  conduction 
band  nonparabolicity  into  the  hydrodynamic  model.  However,  for  this  initial  investigation  with  the  new 
multi-valley  model,  which  incorporates  a  separate  set  of  hydrodynamic  equations  for  each  conduction 
band  (r,  L  and  A)  and  allows  for  electron  exchange  through  energy  dependent  relaxation  parameters, 
the  nonparabolic  correction  terms  will  be  suppressed.  This  was  done  to  simplify  the  analysis  and  to 
provide  a  base  for  future  aoaparabolic  studies.  Additional  features  of  the  model  include  self-consistent 
potentials  aad  realistic  boundary  conditions. 
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Abstract.  In  this  paper,  the  results  of  studying  a  submicron  GaAs  ballistic  diode, 
using  a  new  multi-valley  (T.  L  and  X)  non-parabolic  hydrodynamic  transport 
model,  are  presented.  Numerical  simulations  indicate  that  accurately  including  the 
effects  of  non-parabolicity  in  the  streaming  terms  and  A-space  transfer  in  the 
velocity  and  energy  equations  is  very  important  in  correctly  determining  the 
conductance  of  the  device.  The  existence  and  amount  of  negative  differential 
conductance  was  determined  to  be  strongly  influenced  by  both  of  these  physical 
factors.  Furthermore,  the  sensitivity  of  device  conductance  to  changes  in  the 
thermal  conductivity  is  diminished  significantly  when  non-parabolicity  is 
accurately  incorporated. 


1.  Introduction 

Advanced  macroscopic  transport  models,  based  upon 
moments  of  the  Boltzmann  transport  equation  (an)  as 
opposed  to  a  full  microscopic  description  (e.g,  the  Monte 
Carlo  (mc)  approach),  are  valuable  tools  to  efficiently 
investigate  electron  dynamics  in  submicron  semicon¬ 
ductor  structures  [1].  However,  these  models  must  ac¬ 
curately  reflect  non-equilibrium  and  non-stationary 
electron  physics  before  they  can  be  applied  successfully  to 
a  wide  variety  of  device  operating  conditions.  Presently, 
the  majority  of  macroscopic  electron  device  models  treat 
the  multi-valley  problem  indirectly  (effective  single¬ 
electron  gas  approximation)  and  ignore  the  full  effects  of 
band  non-parabolicity  [2].  When  more  complete  macro¬ 
scopic  models  are  implemented,  they  sometimes  predict 
questionable  results.  For  example,  a  two-valley  (parabo¬ 
lic)  hydrodynamic  model  has  been  used  to  derive 
negative  differential  conductance  in  a  GaAs  fet  while 
corresponding  mc  simulations  indicate  no  such  effect  [3]. 

In  this  paper,  a  new  hydrodynamic  transport  model 
(iitm)  suitable  for  studying  hot  electron  transport  in 
multiple  non-parabolic  conduction  bands  is  investigated. 
The  impact  of  Jc-space  (intervalley)  transfer  and  band 
non-parabolicity  is  examined  by  studying  variations  of 
this  htm.  To  determine  the  influence  of  these  effects 
individually,  these  various  forms  of  the  model  were 
applied  to  a  submicron  GaAs  ballistic  diode.  Numerical 
simulation  results  indicate  that  including  the  effects  of 
non-parabolicity  in  the  streaming  terms  and  intervalley 
scattering  in  the  velocity  and  energy  equations  strongly 
influences  the  characteristics  of  the  device  in  the  non¬ 
linear  (saturation)  region.  These  results  demonstrate  that 


correctly  including  non-parabolicity  effects  in  the  hydro- 
dynamic  model  reduces  the  amount  of  negative  dif¬ 
ferential  resistance  (ndr).  Also,  the  amount  of  nor  was 
found  to  be  less  sensitive  to  changes  in  the  thermal 
conductivity  when  non-parabolicity  is  accurately  in¬ 
corporated.  Therefore,  these  results  verify  the  importance 
of  both  non-parabolicity  and  intervalley  effects  in  the 

HTM. 


2.  The  hydrodynamic  transport  model 

For  these  studies,  a  submicron  GaAs  ballistic  diode  (see 
inset  in  figure  I)  will  be  considered  with  transport  in  each 
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Figure  1.  Total  electron  currant  density  against  applied 
electric  potential  for  a  submicron  GaAs  ballistic  diode  (see 
inset).  Results  ere  for  variations  in  the  htm  (broken 
curves)  end  s  mc  study  (full  curve).  Here  1*1  in  the 

HTM. 
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ABSTRACT 


A  novel  real-apaca  transfer  logic  tranaietor  (RSTLT)  b 
propoaad  baaed  on  the  real- apace  traatfer  of  hot  clectrona  in 
a  four-terminal  hetcrojeactioo  mkrostructurc.  Sdf-coubt- 
ent  steady-state  and  tranaieat  Monte  Carlo  rimnlatiena  dam* 
onatrate  that  the  RSTLT  features  both  logic  flexibility  and 
ultra- last  switching  speed.  Calculated  results  show  that  the 
proposed  RSTLT  realises  NOT  and  EQUIVALENT  logic 
functions  in  a  single  hetcro junction  device  and  a  conser¬ 
vative  estimate  of  the  characteristic  switching  tune  is  ~  3 
pace. 


INTRODUCTION 


Real-space  transfer  (RST)  [1]  b  the  thermionic  emis¬ 
sion  of  hot  electrons  (from  one  semiconductor  layer)  over 
hctcroiaterfrce  barriers  to  different  semiconductor  layers. 
Successful  implementations  of  the  RST  concept  have  been 
demonstrated  by  experimental  realisations  of  the  negative 
resistance  field-effect  transistor  (NERFET)  [2],  charge  in¬ 
jection  transistor  (CHINT)  (3)  in  both  AlGaAs/GaAs  and 
InGaAs/InAlAs  material  systems,  and  real- apace  transfer 
transistor  (RSTT)  (4)  using  strained  InGaAs/AJGaAs/GaAs 
layers.  Promising  performances  have  been  achieved  in  these 
devices,  with  reported  values  of  traasconductancs  exceed¬ 
ing  2300  mS/mm  (3)  and  unity-current  gat?  frequency  of 
60  GHs  |4)  at  room  temperature.  Abo,  a  NORAND  func¬ 
tional  element  baaed  on  a  multi-terminal  CHINT  structure 
has  been  proposed  recently  [5],  demonstrating  the  poten¬ 
tial  logic  applications  of  real-space  transfer  devices.  In  this 
work,  we  propose  a  novel  logic  circuit  element  •  the  real- 
space  transfer  logic  transistor  (RSTLT).  The  operation  of 
this  proposed  mkrostructure  b  based  on  the  concept  of  hot 
electron  real-space  transfer  and  the  fret  that  the  spatial 
location  of  electron  RST  b  determined  by  applied  bias  and 
hcteroiaterface  energy  barrier  height.  Sefrconsbtcat  ensem¬ 
ble  Monte  Carlo  simulations  are  used  to  theoretically  'build’ 
a  RSTLT  with  flexible  logic  functions  and  ultra-bat  speed. 


RSTLT  STRUCTURE  AND  SIMULATION 
METHOD 


Figure  1  shews  the  proposed  RSTLT  structure.  The  de¬ 
vice  has  two  separate  collectors  (Cl  and  C2)  as  output  ter¬ 
minals  which  are  connected  through  load  dements  (R)  to 
the  positive  power  supply  (Vfc),  a  drain  (D)  input  termi¬ 
nal  and  a  common  source  (S)  terminal  which  b  at  ground 
potential.  When  D  b  biased  at  logic  low  (positive  logic), 
channel  electrons  confined  in  the  potential  wdl  between  the 
channel- barrier  heterointerfrees  will  not  have  enough  energy 
to  overcome  the  energy  barrier  until  reaching  C2  where  they 
are  collected.  In  this  case  Cl  b  logic  frph  and  C2  b  logic 
lew  (terminal  voltages  of  Cl  and  C2  are  high  and  low,  re¬ 
spectively).  On  the  other  hand,  if  D  b  high,  most  electrons 
gain  sufficient  energy  to  overcome  the  barrier  while  pessing 
under  Cl.  Therefore,  Cl  becomes  low  and  C2  hsph.  Thus, 
the  RSTLT  realises  the  logic  NOT  (Cl)  and  EQUIVALENT 
(C2)  functions  ef  the  input  contra!  signal  (D),  Le.,  Cl  ■  77, 
C2  «  D,  as  shewn  in  the  inset  of  Fig.  1.  In  order  to  ensue 


Figure  1.  Proposed  RSTLT  structure.  The  inset  on  the  top 
right  shows  the  logic  symbol  of  the  structure. 


CHJ975-9'91/0000-0515  S1.00  C  1991 IEEE 


18.6.1 


IEDM  91-515 


of  Seaiconduccor  Devices  and  Integrated  Circuits  (Hay,  199Z,  Vienna,  Austria; 


AN  EXTENDED  SCHARFETTBR-CUMMEL  ALGORITHM  SUITABLE 
FOR  SOLVING  A  NONPARABOLIC  HYDRODYNAMIC  TRANSPORT  MODEL 

D.L.  WOOLARD,  T.A.  WINSLOW,  M.A.  LITTLEJOHN  and  RJ.  TREW 

Electrical  and  Computer  Engineering  Department 
North  Carolina  State  University 
Raleigh,  North  Carotiaa  2TSW-7tIt  U.S.A. 


Abetract  •  Thie  paper  prevent*  an  extended  Schaifetter-Gummel 
carrent-deneity  discretisation  forma  la  tuIBcient  for  *oleing  an  ad¬ 
vanced  hydrodynamic  electron  traneport  model.  Thb  Improved  die- 
cretixatioa  eapreeshm  la  specifically  developed  far  a  new  hydrody¬ 
namic  traneport  model  suitabk  far  rtudytag  hot-electron  traneport 
in  vemicond acting  materials  with  nonparahellc  condactlon  band*. 
While  this  attended  forma  lotion  properly  incorporate*  many  ad¬ 
vanced  effect*  associated  with  hot -electron  transport,  its  primary 
/rater*  I*  the  accurate  treatment  o f  the  ran  vet  ilea  term  (or  displace¬ 
ment  energy  term)  which  Is  present  in  hath  conventional  (La.  drifted 
MaxwcUtou  model)  and  the  advanced  nonpar abofle  model  dlstonid 
In  this  paper.  Oor  derivation  overcomes  nonlinear  problem*  intro¬ 
duced  by  this  convection  term  and  resalts  in  an  Improved  expression 

to  approximate  current  density. 


I.  INTRODUCTION 


At  prevent,  Monte  Carlo  (MC)  mathode  and  hydrodynamic  traneport 
models  (IITM’s)  are  the  two  mast  popnlar  approaches  utilised  to  in¬ 
vestigate  hot-electron  phenomenon  In  suhmicroa  structures.  While 
MC  methods)!),  which  solve  the  full  Boitimana  transport  agnation 
(DTE),  can  easily  Incorporate  detailed  traneport  physics,  thb  ap¬ 
proach  can  become  Impractical  in  some  instances  doe  to  the  extensive 
computational  requirements.  RTM'i,  which  are  based  upon  macro 
acopic  conservation  equations  derived  from  m  aments  af  the  8TB,  in¬ 
clude  nonet et ionary  and  hot  electron  effects)])  while  reqnlrlng  much 
lee*  computation  time  to  i*a crate  solutions.  The  physical  and  prac¬ 
tical  attributes  of  HTM'ehti*  land  to  Uteb  extanaiva  application.  In 
turn,  advanced  IfTM'a  have  bean  developed  which  are  accurate  hr 
various  device  structures  and  operating  condition *|3j.  While  IITM’s 
represent  a  considerable  reduction  In  complexity  (a*  compared  to 
solving  the  full  BTE),  their  solution  far  arbitrary  desks  structure* 
Is  still  a  very  formidable  mat  hematics  problem.  When  IITM's  are 
applied  aelf-consietently  with  Posslau**  squ alien,  one  must  solve  a 
set  of  noaHnear  coupled  partial  differential  squat  lone.  Thb  tasfc  re¬ 
quires  the  application  of  accurate  and  (tabic  numerical  techniques. 
Since  the  success  of  any  numerical  stamtotka  b  critically  dependent 
on  convergence  rate  and  numerkel  stability,  the  choice  of  an  efficient 
discretisation  scheme  b  very  Important)!). 

It  b  well  known  that  a  central  difference  approach  b  not  adequate  for 
d  beret  king  the  current  continuity  equation  due  to  the  development 
of  esclllalians  or  “wlyffee**  In  the  point  ion) I).  To  Insure  numerical 
stability,  Scharfstter  and  Oummsfffij  Introduced  a  physics- based  cur- 
rent  density  faemalatlou  be  the  drift  dfffcriou  (PD)  model  (l.e.  they 
ueenmed  current  density  voetas  weakly  with  position).  The  success 
the  Schuefetter-  Ournmul  (SO)  schema.  In  DD  model  elmubtbni,  has 
motivated  various  geueralbatlsns  of  thb  dbcrotbatlea  technique  for 
the  hydrodynamic  equations) 7).  However,  most  ixtvwsluns  of  the 
SG  formula  lor  IITM’s  either  ignore  convention  terms)!)  or  eba  as¬ 
sume  o  priori  that  they  are  doffmd  In  a  previous  step  af  an  Karelian 
procedure)?].  In  thb  paper,  an  analytical  method  b  presented  to  de¬ 
velop  on  extended  SG  current  density  fsrmuk  which  ffiroctly  Indudss 
the  effects  af  the  nouHnsar  convention  term.  Thb  SO  algorithm  b 
spadffcaffy  developed  for  use  In  salving  a  new  uonparabeffc  HTM. 
Abo,  the  extended  fanuab  b  ismparcd  to  one  which  Ignores  energy 
displacement  to  Mustcate  the  bison  reef  thb  neuffnoar  effect  on  the 
predicted  current  density. 


II.  TIIE  NONPARABOLIC  MODEL 


Prevbnely,  a  hydrodynamic  model  suitable  for  studying  bot-viectron 
transport  to  semiconducting  material*  with  nonparabolk  conduction 
bonds  was  presented))).  Thb  new  model  contains  admitted  non- 
parabolic  streaming  (collbion  independent)  terms  which  were  de¬ 
rived  by  applying  a  unique  set  of  moment  operator*  (l.e.  order* 
of  nonlinear  k-*pace  velocity  a*  opposed  to  momentum)  to  the  col- 
Ustonks*  Boltxmana  transport  aquation.  The  resulting  system  of 
moment  (balance)  equation*  were  closed  (terminated)  by  Introduc¬ 
ing  an  empirical  distribution  function  which  accurately  approximate* 
the  properties  of  electron*  drifting  to  non  parabolic  conduction  bond*. 
To  treat  totervally  tranehr  effects  that  are  present  In  multi- vaffey  aye- 
term  (04.  GaAe),  a  throe-valley  (I,  L  and  X )  ensemble  relaxation 
time  modal  ha*  been  developed  and  implemented) 8).  The  full  multi- 
valley  nonparabolk  HTM,  with  a  phsnomcnotogkal  heat  Row  vector 
term,  consist*  of  the  transport  equation* 
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for  each  condactlon  hand  valley  (I  m  T,  L  and  X).  In  the  prevloa* 
model  equations,  q  b  the  electron  concentration,  v/  b  the  average 
eketroa  velocity,  si  b  the  average  electron  energy  and  F  n  fV.d 
b  the  farce  on  the  ebetron.  to  thb  HTM,  the  resulting  conrtlto- 
live  relations  are:  an  average  effective  mas*  **(vj,  *>  m)(t  4 

9oii|T^  -4  T-ltaiiifq)1 4-  cqmf  |v<|*),  an  effective  velocity  prewar* 
|P,|J  »  (m;/rh;)u,*,r^|I).  an  effective  energy  pressure  [P»J  ■ 
*il|?a(l  4-  4-  «tml|vt)*))l|,  a  nonparabaffc  boat  Row  vet- 

lor  qi  <u  •ttl(m|/**)*(fclTri)l(l-i-3lw<(m)/to))Aprwi)n|V<,and  th* 
nonparabolk  energy  partition  relation  aq  a  (3/2)hpT,//(vifTw)  4 
(*;(v,,T^)/I)|v,|*  with  an  effective  temperature  of  T.//(v,,T^) » 
T«r(l  4-  iaikpTmi  4-  ojmf|vi)’).  These  previous  rehflom,  which  ex¬ 
press  the  corrected  streaming  terms  as  functions  of  v,  w  and  new- 
parahoBc  conduction-bond  factor  o  were  verified  far  both  ststlmory 
and  nonet stlonsry  transport  to  GnA*  by  two  ludspsadmt  Meat* 
Carlo  stomlation  *tudies|3).  Features  of  thb  model  which  ffiffer  heat 
th*  conventional  HTM  Include?  (1)  a  non  roust  snt  average  effeetke 
mas*  dl)(vi,To<),  (3)  a  distinctly  different  velocity  and  energy  pres¬ 
sure  (L*.  (Fell  y<  (Fw(|),  and  (3)  a  first- principle*  nansma  tost 
Row  vector  qi.  In  tide  model,  q*  represent*  the  heat  Row  resaltbg 
from  th*  symmetric  campon  ml  of  tbs  velocity  space  electron  dhtri- 
butian  In  ueuparubolle  boards.  ^Therefore,  the  Wlede^nann—  ^bons  pb^^ 
naminslBgltal  heat  Row  term,  qf  to  Eq.  (!),  b  Introduced  to  astsoal 
far  passible  asymmetry  of  the  ebetron  dbtrfbutteu.  to  thb  boat  to 
term,  pn  b  the  near  equilibrium  nwbiDty  and  Q  b  a  muMpRtaHw 
constant  associated  wkh  tb*  particular  material  ustoercaasidsratba. 
Dtasipatlan  and  Intarvalley  transfer  effects  (l.e.  lust  terms  in  Bqa  (I) 

-  (3))arotodudsd  using  an  ensembb  rotor atbn  time  model  Hamid 
to  farm  to  the  one  pmpoeod  by  Blateh)aer)R{.  Te  templets  themedd, 
the  various  energy-dependent  wtoxntlon  time*  fare  wudU-mBeyrys 
tarn  (OnAs)  have  bam  ebtslaad  ham  stationary  MC  cabufattouH 
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Abstract  •  The  effects  of  laterally  varying  doping  profile  on 
device  electron  transport  properties  are  studied  using 
self-consistent  ensemble  Monte  Carlo  simulations.  Significant 
performance  improvement  is  observed  for  N*  —  N  —  N*  struc¬ 
tures  with  ramp- down  doping  and  multi-spike  doping  profiles 
in  the  active  N  region.  The  prospects  for  application  of  later¬ 
ally  varying  doping  schemes  in  field-effect  transistor  structures 
are  discussed. 


The  utilisation  of  'bandgap  engineering’  has  led  to  a  va¬ 
riety  of  novel  improved  device  structures  in  the  past  decade. 
Successful  integration  of  'bandgap  engineering'  with  'doping 
engineering'  provides  additional  flexibility  and  new  methods 
for  the  future  development  of  advanced  semiconductor  devices. 
Currently,  most  semiconductor  devices  employ  techniques  that 
control  doping  profiles  in  the  vertical  dimension  (perpendic¬ 
ular  to  the  direction  of  carrier  transport  using  04.,  uniform 
doping,  modulation  doping,  and  delta-doping).  The  rapid  de¬ 
velopment  of  new  technologies  such  as  focused  ion  beams  (FIB) 
makes  possible  mask-less  direct  ‘writing*  of  controlled  dopants 
into  semiconductors  [1|.  Thus,  the  doping  distribution  can  he 
tailored  in  the  lateral  dimension  (parallel  to  the  direction  of 
carrier  transport).  This  added  degree  of  freedom  for  'doping 
engineering'  can  possibly  result  in  novel  device  structures  with 
improved  performance.  The  purpose  of  this  paper  is  to  present 
results  of  a  theoretical  study  of  the  effects  of  laterally  varying 
doping  profiles  on  device  electron  transport  properties,  which 
supports  this  proposition. 

In  this  work,  self-consistent  ensemble  Monte  Carlo  simula¬ 
tions  are  performed  to  investigate  the  effects  of  laterally  vary¬ 
ing  doping  on  electron  transport.  The  Monte  Carlo  model  em¬ 
ploys  a  three-valley  (r-L-X),  aoa parabolic  band  structure  cou¬ 
pled  with  a  solver  for  Poisson's  equation.  Trajectories  of  a  large 
number  of  sample  electrons  are  traced  in  real  space  and  in  mo¬ 
mentum  space.  Relevant  scattering  mechanisms  (polar  optical, 
equivalent  and  nonequivalent  intervalley,  ionised  impurity,  and 
electron-electron)  are  included.  In  addition,  the  effect  of  degen¬ 
eracy  for  r- valley  electrons  is  taken  into  account  by  employing 
an  approximate  Fermi-Dirac  distribution  with  calculated  local 
quasi-Fermi  level  and  electron  temperature  |2|.  Charge  neu¬ 
trality  in  a  portion  adjacent  to  the  electrodes  is  maintained 
throughout  the  simulation,  which  serves  as  the  criterion  tor 
electron  injection.  Material  and  energy  band  parameters  for 
GaAs  are  the  same  as  those  in  Ref.  |3|.  Lattice  temperature 
is  assumed  to  be  300  K.  More  detailed  model  description  can 
be  found  elsewhere  |4|.  This  transport  study  assesses  advan¬ 
tages,  differences,  and  potential  applications  of  different  later¬ 
ally  varying  doping  profiles  on  device  performance.  The  test 
device  is  a  GaAs  —  ff(OJpm)  —  (0.5pm)  struc¬ 

ture  with  varied  doping  profile  along  the  direction  of  electron 
transport  in  the  N  region,  la  particular,  two  groups  of  laterally 


varying  doping  schemes  are  studied.  They  are:  1)  ramp-doping 
|i-e-,  linearly  increasing  (ramp-up)  or  decreasing  (ramp-down) 
doping  density  within  the  N  region),  and  2)  spike-doping  (i.c., 
introducing  one  or  more  N*  spikes  in  the  N  region).  In  order 
to  have  meaningful  comparisons  among  structures  which  em¬ 
ploy  different  doping  profiles,  we  apply  conditions  of  identical 
minimum  doping  density  and  total  integrated  dopant  in  the  N 
region. 

Any  successful  implementation  of  a  lateral  doping  scheme 
should  result  in  improved  current  drive  capability.  Among  all 
the  devices  studied,  our  simulation  results  reveal  that  the  most 
significant  improvement  can  be  achieved  by  employing  a  struc¬ 
ture  with  a  ramp-down  doping  scheme.  For  the  same  applied 
bias,  the  ramp-down  doped  structure  exhibits  more  than  50  % 
overall  improvement  in  current  drive  capability  compared  to 
that  for  the  structure  with  the  ramp-up  doping  scheme.  This 
is  shown  in  Fig.  1,  where  current-voltage  characteristics  of 
the  structures  with  ramp-up,  ramp-down,  and  constant  mini¬ 
mum  N  doping  profiles  are  given.  The  conduction  band  profile 
(increasing  exponentially  along  the  direction  of  electron  trans¬ 
port)  generated  by  the  ramp-down  doping  scheme  and  the  re¬ 
sultant  accelerating  electric  field  profile  provide  favorable  dec- 
iron  transport  conditions  in  which  electrons  accelerate  more 
efficiently.  Fligh  average  electron  velocity,  low  average  electron 
energy,  and  high  r- valley  electron  occupancy  are  obtained  for 
the  ramp-down  doped  structure.  These  results  indicate  that, 
on  average,  mote  electrons  in  the  ramp-down  doped  structure 
accelerate  while  they  stay  in  the  T  valley.  Results  from  spike- 
doped  structures  show  that  introducing  an  N*  spike  in  Ike 
active  electron  transport  (N)  region  brings  improved  device 


Figure  I.  Current- voltage  characteristics  for  the  ramp-doped 
N*  -  N  —  N*  structures.  Doping  density  varies  Bneariy  be¬ 
tween  5  x  10'*em-*  and  5  x  10,Tem-*.  Result  for  N*  -  N - 
N*  structure  with  constant  (minimum)  N  doping  density  of 
5  x  I0'*cm~a  is  induded  for  comparison. 
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1.  Introduction; 

Enhancement  in  the  high  frequency  performance  of  transistors  is  dependent 
both  on  improved  materials  properties  and  smaller  transit  distances.  This  has 
resulted  in  improved  performance  as  devices  have  progressed  from  MESFETs 
(Wang  etal)  to  GaAs  based  pseudomorphic  HEMTs(Tan  etal)  to  InP  based 
lattice-matched (Mishra  etai,  Chao  etal)  and  pseudomorphic(Thompson  etal) 
HEMTs.  The  GalnAs  based  devices  have  demonstrated  foe  highest  speeds 
because  of  a  combination  of  higti  electron  mobility,  peak  velocity  and  sheet 
charge  density  in  the  GalnAs  channel  of  an  AllnAs-GalnAs  HEMT.  The 
reduction  in  gate  length  to  50  run  has  increased  foe  extrinsic  fr  of  the  HEMT 
to  292  GHz  at  room  temperature(Thompson  etal)  but  the  fmax  of  the  HEMT 
has  been  restricted  by  foe  rapid  increase  in  the  output  conductance  at  the 
small  gate  length.  The  output  conductance  is  caused  by  a  combination  of 
substrate  injection  and  channel  length  modulation.  Attempts  to  reduce  foe 
substrate  injection  by  using  a  combination  of  p-doped  and  wide  band  gap 
buffer  layers  have  had  limited  success.  In  this  study,  we  investigate  the  effects 
of  doping  foe  channel  of  an  AllnAs-GalnAs  modulation  doped  transistor  n 
and  p  type  and  evaluate  foe  effect  cm  the  DC  and  RF  characteristics  of  devices 
with  0.25  pm  gate  length.  Devices  with  1  pm  gate  length  were  also  studied  to 
determine  foe  gate  length  dependence.  The  motivation  is  to  evaluate  foe 
effect  of  foe  two  doping  types  on  the  distribution  of  electric  held  in  the 
channel  and  study  its  effect  on  output  conductance  and  electron  transport. 
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Figure  1:  Schematic  of  the  Layer  Structures  Used  in  foe  Study. 
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Two-Dimensional  Analysis  of  Short-Channel 
Delta-Doped  GaAs  MESFET’s 
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Salah  M.  Bedair,  and  Larry  C.  Witkowski 


Abstract— Key  design  parameters  for  delta-doped  GaAs 
MESFET’s,  such  as  delta-doping  profile,  top  layer  background 
doping  density,  and  scaling  of  lateral  feature  size,  are  investi¬ 
gated  using  a  two-dimensional  numerical  simulation.  A  three- 
region  (delta-doped  conducting  channel,  top  layer,  and  sub¬ 
strate)  velocity-field  relation  is  implemented  in  the  model  as 
appropriate  for  tbe  particular  device  structure  which  is  simu¬ 
lated.  Simulation  results  show  excellent  agreement  with  our 
fabricated  0.5- pm  gate-length  delta-doped  GaAs  MESFET’s 
based  on  atomic  layerepitaxy  material.  Anextrinsictranscondisc- 
tance  of  370  mS/mm  and  a  drain-source  current  of  270 
mA/mm  are  obtained  for  typical  devices,  and  the  maximum 
transconductance  is  as  high  as  400  mS/mm,  which  are  the  best 
dc  results  yet  reported  for  0.5-pm  gate-length  delta-doped  GaAs 
MESFET’s.  Considerations  of  design  and  optimization  are  dis¬ 
cussed  in  terms  of  threshold  voltage  sensitivity,  transconduc¬ 
tance,  current  drive  capability,  and  cutoff  frequency,  based  on 
both  simulation  and  experimental  results. 


I.  Introduction 

THE  DRIVE  to  achieve  high-speed  and  high-fre¬ 
quency  electronic  systems  has  led  to  continued  efforts 
to  develop  new  device  structures.  One  veiy  promising 
structure  is  a  field-effect  transistor  (FET)  which  incorpo¬ 
rates  the  delta-doping  technique.  Advanced  delta-doped 
FET  structures  such  as  GaAs  MESFET’s  [1J— £41, 
HEMT’s  on  GaAs  [5]-[7]  and  InP  substrates  [8],  as  well 
as  pseudotnorphic  HEMT’s  [9]  have  been  demonstrated. 
The  advantages  of  delta  doping  in  these  device  structures 
include  high  current  drive  capability,  high  transconduc- 
tance,  improved  threshold  voltage  control,  and  improved 
breakdown  characteristics. 

Delta-doped  GaAs  MESFET’s  based  on  materials 
grown  by  MBE  and  MOCVD  materials  have  been  inves¬ 
tigated  experimentally  by  several  research  groups  [1]-[4J. 
While  the  results  are  promising,  they  also  show  a  wide 
diversity  and  fall  short  of  theoretical  predictions  [1].  The 
reported  maximum  transconductance  of  0.5- pm  gate- 
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length  delta-doped  MESFET’s,  for  example,  ranges  from 
75  mS/mm  [1]  to  140  mS/mm  [2),  to  die  400  mS/mm 
transconductance  reported  in  this  paper.  Schubert  et  al. 
[1]  studied  properties  of  delta  doping  in  GaAs  and  pre¬ 
dicted  that  the  maximum  intrinsic  transconductance  of  a 
delta-doped  MESFET  could  be  higher  than  that  of  a  com¬ 
parable  AlGaAs/GaAs  HEMT  in  the  gate  length  range  of 
0.1  to  2.0  pm.  Differences  between  experimental  results 
and  theoretical  predictions  may  stem  from  unknown  ma¬ 
terial  properties  and  differences  or  uncertainty  in  process¬ 
ing  conditions.  On  the  other  hand,  design  and  optimiza¬ 
tion  of  such  devices  play  an  important  role  in  guiding  the 
improvement  of  device  performance.  It  is  of  practical  im¬ 
portance  to  identify  the  key  parameters  and  the  tradeoffs 
needed  for  optimized  device  performance  in  specific  ap¬ 
plications. 

The  intent  of  this  paper  is  to  present  a  generalized  anal¬ 
ysis  of  delta-doped  submicrometer  MESFET’s.  We  use  a 
two-dimensional  drift-diffusion  model  to  simulate  the  in¬ 
fluence  of  key  design  parameters  on  delta-doped  MES¬ 
FET  performance.  Experimental  0.5-nm  gate-length 
delta-doped  MESFET’s  grown  by  atomic  layer  epitaxy 
(ALE)  have  also  been  fabricated.  Excellent  agreement  be¬ 
tween  simulation  and  experimental  measurement  has  been 
obtained.  An  extrinsic  transconductance  of  370  mS/mm 
and  a  drain-source  current  of  270  mA/mm  are  measured 
for  typical  devices,  while  the  maximum  transconductance 
is  as  high  as  400  mS/mm.  These  are  the  best  experimen¬ 
tal  results  yet  reported  for  0.5-/tm  gate-length  delta-doped 
GaAs  MESFET’s.  Based  on  the  agreement  between  sim¬ 
ulation  results  and  experimental  results  for  0.5-pm  de¬ 
vices,  we  then  analyze  performance  dependence  on  key 
device  parameters  and  outline  some  design  considerations 
and  tradeoffs  for  achieving  improved  delta-doped  MES¬ 
FET  operation. 

II.  Device  Fabrication 

A  schematic  illustration  of  a  delta-doped  GaAs  MES¬ 
FET  is  shown  in  Fig.  1(a).  A  very  narrow  doping  profile 
(Fig.  1(b))  is  sandwiched  between  two  undoped  (or  rela¬ 
tively  low-doped)  layers.  The  device  structure  and  the 
doping  profile  in  the  figures  are  representative  of  those 
employed  in  the  experimental  device,  although  the  actual 
device  fabricated  for  measurement  in  this  paper  contains 
a  recessed  gate.  Tbe  doping  profile  used  in  the  simulation 
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Abstract— Ensemble  Monte  Carlo  simulations  are  employed 
in  order  to  explore  the  feasibility  of  a  novel  real-space  transfer 
logic  transistor  (RSTLT)  structure.  The  operational  principles 
of  the  proposed  RSTLT  are  based  on  the  concept  of  hot  electron 
real-space  transfer  (RST),  including  the  fact  that  the  spatial 
location  of  electron  RST  is  determined  by  applied  bias  and  het¬ 
erointerface  energy  barrier  height  in  a  multiterminal  hetero- 
junction  microstructure.  The  results  of  two-dimensional,  self- 
consistent  steady-state  and  transient  simulations  demonstrate 
that  the  proposed  RSTLT  features  ultrafast  current  switching 
which  can  be  used  to  realize  NOT/EQUIVALENT  logic  func¬ 
tions  in  a  single  beterojunction  device.  The  logic  operation  is 
easily  extended  to  NOR/AND  functions.  A  conservative  esti¬ 
mate  of  the  characteristic  delay  time  for  current  switching  is 
-3  ps  in  the  proposed  RSTLT  structure. 


I.  Introduction 

THERMIONIC  emission  of  hot  electrons  (from  one 
semiconductor  layer)  over  heterointerface  banier(s) 
to  different  semiconductor  layers  in  the  presence  of  high 
electric  fields  is  termed  real-space  transfer  (RST)  [1]. 
Successful  implementations  of  the  RST  concept  have  been 
demonstrated  by  experimental  realizations  of  the  nega¬ 
tive  resistance  field-effect  transistor  (NERFET)  [2]-[4], 
the  charge  injection  transistor  (CHINT)  in  both  Al- 
GaAs/GaAs  [5],  [6]  and  InGaAs/InAlAs  [7],  [8]  mate¬ 
rial  systems,  and  the  real-space  transfer  transistor  (RSTT) 
using  strained  InGaAs  / AlGa As  / GaAs  [9]  heterostruc¬ 
tures.  Promising  performances  in  both  dc  and  microwave 
applications  have  been  achieved  in  these  devices,  with  re¬ 
ported  transconductance  values  exceeding  2300  mS/mm 
[6]  and  unity-current  gain  frequency  of  60  GHz  [9]  at 
room  temperature.  Recently,  a  CHINT  logic  element  to 
realize  nor/ and  functions  [10]  has  been  proposed,  dem¬ 
onstrating  the  potential  logic  applications  of  real-space 
transfer  devices.  Theoretical  investigations  of  electron 
transport  properties  and  related  device  physics  in  real- 
space  transfer  structures  have  also  been  performed  [11]— 
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[18],  These  studies  provide  a  foundation  for  device  struc¬ 
tures  with  novel  operating  principles. 

It  is  worth  noting  that  in  a  real-space  transfer  device 
(such  as  the  CHINT),  hot  electron  RST  is  a  localized  pro¬ 
cess  in  which  the  spatial  occurrence  of  RST  is  determined 
by  factors  such  as  the  device  configuration,  heterojunc¬ 
tion  material  parameters,  and  operating  bias  conditions. 
This  phenomenon  has  been  discussed  by  interpreting  ex¬ 
perimental  I-V  characteristics  of  CHINT  structures  using 
positive  feedback  initiated  by  a  local  ‘hot  spot’  in  the  de¬ 
vice  channel  [8].  Also,  Monte  Carlo  analysis  of  three- 
terminal  RST  devices  has  indicated  the  localized  nature 
of  electron  RST  [13].  Recently,  we  have  proposed  a  novel 
real-space  transfer  logic  transistor  (RSTLT)  [19]  based  on 
ensemble  Monte  Carlo  studies  of  hot  electron  transport  in 
a  four-terminal  heterojunction  microstructure.  The  oper¬ 
ation  of  the  proposed  RSTLT  is  based  on  the  fact  that  the 
spatial  location  of  electron  RST  is  determined  primarily 
by  applied  bias  and  heterointerface  energy  barrier  height. 
The  dc  and  transient  current  switching  characteristics  of 
the  proposed  RSTLT  have  been  discussed  briefly  [19] 
(here  we  define  current  switching  as  the  current  exchange 
among  device  terminals  while  transient  switching  times 
for  the  intrinsic  device  are  denoted  as  characteristic  delay 
times).  We  have  demonstrated  the  the  potential  logic  ap¬ 
plications  based  on  multiple  terminal  current  switching 
characteristics  of  the  RSTLT. 

In  this  work,  detailed  theoretical  investigations  using 
self-consistent,  two-dimensional  ensemble  Monte  Carlo 
simulations  are  presented  which  focus  on  exploring  the 
feasibility  of  the  RSTLT  structure  in  digital  applications. 
The  proposed  structure  differs  from  the  CHINT  logic  ele¬ 
ment  [10],  i.e.,  the  RSTLT  employs  a  different  device 
configuration  and  terminal  arrangement  to  achieve  shorter 
transient  delay  time  while  maintaining  logic  flexibility. 
This  provides  another  attractive  approach  to  realize  ultra¬ 
fast  multilogic  functions  in  a  single  heterojunction  de¬ 
vice.  This  computer  experimental  study  is  used  to  theo¬ 
retically  “build"  a  device  which  allows  the  logic 
flexibility  and  fast  speed  of  available  real-space  transfer 
devices  to  be  explored.  Electrical  characteristics  and  re¬ 
lated  device  physics  of  the  proposed  RSTLT  structures 
are  examined  in  terms  of  material  parameters,  device  fea¬ 
ture  size,  and  bias  conditions.  In  the  next  section,  we  de- 
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Theoretical  results  of  electron  transport  in  n-type  In0,52A]04(As  are  presented.  The  transport 
properties  of  this  important  semiconductor  were  obtained  using  the  Monte  Carlo  method.  In 
particular,  velocity-electric  field  characteristics  for  different  temperatures  and  doping 
concentrations  in  bulk  Ina}2Al0.4|As  are  calculated  for  the  first  time.  Physical  parameters  for 
In^Al^As  (which  is  lattice-matched  to  InP  and  In0]]Ga0.47As)  were  obtained  based  on 
interpolation  of  available  experimental  and  theoretical  results  for  InAs,  A  LAs,  and 
In,,  75AJ,)  uAs.  Our  study  suggests  that  Ioq  S2AIo.4«As  has  electron  transport  properties  which  are 
comparable  to  and  complimentary  with  those  of  other  materials  lattice-matched  to  InP. 


Heterojunction  devices  based  on  the  lattice-matched 
Ina32A](U(As/ln0.5}Ga0l47As/lnP  material  system  have 
been  the  focus  of  extensive  experimental  studies  due  to 
their  promising  potential  for  digital,  microwave,  and  opti¬ 
cal  applications.  Impressive  device  performance  has  been 
achieved  for  the  high  electron  mobility  transistor 
(HEMT),1  charge  injection  transistor  (CHINT),2  and 
semiconductor  laser*  fabricated  from  this  material  system. 
However,  certain  material-related  issues  remain  which 
merit  careful  study  if  potential  applications  of  these  devices 
are  to  be  fully  exploited.  As  an  example,  high  output  con¬ 
ductance  has  been  repeatedly  observed  in  deep-submicron 
( £0.3  /tm)  channel-length  HEMTs  using  the 
In0  ]2Alo  4|As/In0  j]Gao  47As/InP  material  system. 1,4,5  This 
performance  degradation  has  been  suggested  to  be  the  re¬ 
sult  of:  ( 1 )  electron  conduction  in  the  In^AI^As  buffer 
layer,  (2)  deep  level  traps  in  the  top  Ino^AlouAs  donor 
supply  layer,  and  (3)  weak  impact  ionization  in  the 
Ino  5jGao  47As  channel  layer.  Studies  on  the  CHINT  have 
also  shown  that  electron  transit  time  in  the  In^Al^As 
collector  barrier  is  the  main  factor  that  determines  the 
total  device  intrinsic  delay  time.2  These  observations  indi¬ 
cate  the  importance  of  achieving  a  better  understanding  of 
electron  transport  in  In^Al^As.  Unfortunately,  there  is 
very  little  available  information  on  material  parameters  or 
electron  transport  properties  for  Iiiq  j2A104IAs.  This  is  a 
noticeable  void  since  there  have  been  extensive  experimen¬ 
tal  and  theoretical  studies  of  InP  and  Ino  jjGa^As,  and 
data  on  In0)2Al04(As  are  required  in  order  to  perform 
thorough  and  more  rigorous  device  simulations  in  the 
InQ  ]2AI04|As/Ino.})Ga047As/lnP  material  system. 

In  this  letter,  we  employ  a  Monte  Carlo  model  to  study 
electron  velocity-electric  field  characteristics  in  n-type 
In032Al04(As  for  different  temperatures  and  doping  con¬ 
centrations.  The  set  of  physical  parameters  for 
Ino32Alo.4|As  required  for  this  study  has  been  obtained 
based  on  linear  interpolation  of  available  data  for  the  InAs 
and  AlAs  binary  components,  while  using  the  few  available 
(theoretical)  material  properties  of  Ina73AI023As  as  a 
guide.  The  material  parameters  which  we  provide  for 
In0  32Al04JAs  and  the  results  of  the  fundamental  electron 
transport  study  will  serve  to  promote  further  experimenta¬ 
tion  and  characterization  of  I no  uAl^ As-related  devices. 


The  Monte  Carlo  method  used  in  this  work  comprises 
a  model  for  the  motion  of  sample  electrons  in  a  three  valley 
(T-L-JT)  semiconductor  with  nonparabolic  analytical  band 
structure  under  the  influence  of  a  uniform  applied  electric 
field.6,7  Relevant  scattering  mechanisms  are  polar  optical 
phonon  scattering,  equivalent  and  nonequivalent  interval¬ 
ley  scattering,  ionized  impurity  scattering,  alloy  scattering, 
and  impact  ionization.  In  the  present  work,  a  single 
longitudinal-optical  (LO)  phonon  mode  is  used.  There  is  a 
possibility  that  two  LO  phonon  modes  due  to  the  InAs-like 
and  AlAs-like  vibrations  could  be  active  in  In032Al04|As, 
since  it  is  known  that  this  condition  does  exist  in  some 
ternary  materials.  However,  this  is  a  subtle  point,  and  de¬ 
tailed  properties  of  the  phonon  spectra  in  this  material  are 
not  sufficiently  known  at  present  to  justify  the  use  of  two 
LO  phonons  in  the  model.  The  use  of  a  single  LO  phonon 
with  an  energy  of  39.5  meV,  which  was  interpolated  from 
local  mode  energies  of  29.6  meV  for  InAs  and  50.1  meV  for 
A  LAs,  should  not  detract  from  the  main  conclusions  of  this 
work. 

The  ionized  impurity  scattering  rate  used  in  our  sim¬ 
ulations  is  based  on  Ridley’s  statistical  scattering  model.* 
In  this  model,  the  scattering  rate  in  terms  of  electron  wave 
vector  k  is  given  as 


1  gy(fc>[.  /  d  ) 

Tmp(k)  d  r,(k)rB„(k)j]’ 


(1) 


where  v((k)  is  the  electron  group  velocity  and 
d[  =  (2vNdop)~xn]  is  the  average  distance  between  impu¬ 
rities.  The  electron  density  is  assumed  to  be  equal  to  the 
doping  concentration  (A^).  The  scattering  rate  l/rBH, 
which  is  derived  from  the  Brooks-Herring  model,9  was 
calculated  by  following  Rucb  and  Fawcett.10 

The  alloy  scattering  rate  is  obtained  using6 
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dy(c)  , 
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where  Al/  is  the  alloy  scattering  potential,  x  is  the  alloy 
composition,  fl  is  the  primitive  cell  volume,  and  y(t) 
=f  ( 1  +at)  for  the  nonparabolic  conduction  band.  We  use 
the  electron  affinities  of  AlAs  and  InAs  to  predict  a  value 
for  A U.*  Here,  5(a)  is  an  energy-dependent  parameter 
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